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EXECUTIVE  SUMMARY 


This  report  presents  a  full  characterization  of  an  optical  retroreflector  array  designed  for  use 
with  satellite  laser  ranging  on  low  Earth  orbit  (LEO)  spacecraft.  The  array  was  designed  to  provide 
a  robust  optical  link  from  spacecraft  in  a  circular  orbit  of  about  1,100  km  using  typical  NASA-like 
ground  stations,  including  the  small  transportable  laser  ranging  systems.  The  array  will  provide 
unambiguous  returns  for  elevation  angles  above  20°  for  daytime  and  nighttime  ranging.  Consistent 
with  requirements  for  precision  position  estimation,  the  optical  configuration  will  provide  phase 
errors  no  greater  than  those  equivalent  to  a  centimeter  in  ranging  errors.  The  nominal  site  location 
for  the  analysis  was  Midway  Research  Center  (MRC)  in  Quantico,  Virginia. 

Based  on  the  analysis  and  experimental  verification  presented  in  this  report,  we  show  that  the 
NRL  LEO  retroreflector  array  will  meet  these  operational  bounds  well  within  margin.  Furthermore, 
the  array  will  also  close  a  link  for  a  system  as  modest  as  the  field-transportable  laser  radar  station 
(FTLRS)  13-cm  aperture  system  on  Capraia  Island  off  the  Tuscan  coast.  This  latter  analysis 
supports  the  potential  of  configuring  a  compact,  transportable  laser  radar  to  obtain  sub-meter, 
near  real-time,  satellite  ephemerides. 

In  the  course  of  the  analysis  and  design,  a  suite  of  tools  was  developed  to  analyze  performance 
for  single  cubes  and  retroreflector  arrays  of  any  size  and  configuration.  The  tools  include  numerical 
models  that  predict  a  given  configuration’s  optical  characteristics  in  terms  of  laser  radar  cross 
section  (LRCS  or  <tlrcs)*  A  second  part  of  the  tool  set  combines  LRCS  with  ground  station 
characteristics,  atmospheric  properties,  and  orbit  dynamics  to  yield  the  numbers  of  both  photons 
and  photoelectrons  for  a  given  orbit  over  a  specific  ground  site. 

The  array  itself  is  compact  and  lightweight;  it  employs  22  1-cm  retroreflectors  mounted  on  a 
hemisphere  to  provide  a  ctlrcs  of  greater  than  104  m2  for  all  elevation  angles  above  20°  with  108° 
field  of  view  (FOV).  It  weighs  221  grams  and  measures  82  mm  in  diameter  by  43  mm  in  height. 

Analysis  was  verified  experimentally  by  using  a  compact  far  field  test  bed.  Comparison  showed 
good  agreement  between  the  array’s  numerically  modeled  and  measured  optical  characteristics. 
Manufacturing  defects  in  the  wave  front  quality  of  the  retroreflectors  direct  energy  to  a  slight 
extent  into  sidelobes  (less  than  2%).  This  redirection  of  energy  into  the  sidelobes  actually  improves 
the  link’s  performance  by  providing  some  compensation  for  velocity  aberration. 

Results  were  compared  to  the  predicted  performance  of  a  single  retroreflector.  The  comparison 
showed  that  the  FOV  would  be  severely  restricted  by  using  the  single  retroreflector,  hence,  orbital 
sampling  would  be  significantly  reduced.  The  number  of  passes  as  a  function  of  elevation  over  a 
LEO  ground  repeat  track  is  shown  and  illustrates  quantitatively  how  limiting  FOV  in  this  manner 
impacts  pass  yield. 

Predicted  performance  in  the  near-infrared  at  1064  nm  is  also  presented.  Advantages  of  this 
wavelength  include  better  transmission  through  the  atmosphere  and  covertness  due  to  transmission 
in  the  nonvisible  region  of  the  spectrum. 

The  array  was  space-qualified  and  testing  specifications  are  given  in  Appendix  C. 
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PERFORMANCE  CHARACTERISTICS  OF  A  RETROREFLECTOR 
ARRAY  OPTIMIZED  FOR  LEO  SPACECRAFT 


1  INTRODUCTION 

This  report  presents  a  full  characterization  of  an  optical  retroreflector  array  designed  for  use 
with  satellite  laser  ranging  (SLR)  on  low  Earth  orbit  (LEO)  spacecraft.  The  array  was  designed  to 
provide  a  robust  optical  link  from  a  spacecraft  in  a  circular  orbit  of  about  1,100  km  using  typical 
NASA-like  ground  stations,  including  a  small  transportable  laser  ranging  system  (TLRS).  The 
optical  configuration  was  to  provide  phase  errors  no  greater  than  those  equivalent  to  a  centimeter 
in  ranging  errors.  The  nominal  site  location  for  the  analysis  was  Midway  Research  Center  (MRC) 
in  Quantico,  Virginia. 

Satellite  laser  ranging  provides  a  powerful  data  type  for  precise  orbit  determination.  Position 
estimation  based  on  direct  detection  SLR  can  have  performance  comparable  to  a  differential  global 
positioning  system  (GPS)  estimation  and  is  used  as  the  referenced  “truth”  by  the  scientific  com¬ 
munity  for  geoscience  and  navigation.  Although  SLR  data  can  generate  a  highly  precise  orbit 
estimator,  it  is  weather-dependent.  Therefore,  this  data  type  is  uniquely  suited  for  independent 
system  performance  validation  of  onboard  spacecraft  systems  and  for  periodic  calibration  [1]. 

Significant  information  given  in  this  report  includes: 

1.  Link  analyses  providing  the  expected  on-orbit  performance  of  the  NRL  LEO  array  for  two 
different  ground  site  telescope  configurations  at  MRC; 

2.  Link  analyses  of  a  single  retroreflector  and  the  NRL  LEO  array  for  a  TOPEX/Poseidon  pass 
over  the  Tuscan  island  of  Capraia  using  a  13-cm  aperture  SLR  system; 

3.  Supporting  analyses  including  the  direct  numerical  computation  of  a  single  retroreflector  and 
the  NRL  LEO  array’s  optical  properties; 

4.  Pass  yield  as  a  function  of  elevation  angle  and  the  impact  of  field  of  view  (FOV); 

5.  Laboratory  measurements  of  the  array’s  optical  properties  that  validate  numerical  analyses; 
and 

6.  Space  qualification  results  of  the  array’s  mechanical  properties. 

The  report  opens  with  a  review  of  SLR  systems  establishing  general  terminology.  As  part  of  the 
review,  the  specific  assumptions  relevant  to  the  link  analyses  are  stated  in  Section  2.2.  Section  3 
describes  the  basic  mechanical  properties  of  the  array.  Section  4  covers  the  numerical  computations 
required  in  the  link  analyses;  notably,  Section  4.1  considers  single  retroreflector  optical  properties 
and  Section  4.2  extends  the  methods  for  a  retroreflector  array.  The  core  results  of  this  report  are 
the  on-orbit  performance  predictions  in  Section  5.  In  Section  6,  the  report  presents  experimental 
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Fig.  1  -  Satellite  laser  ranging  (SLR)  illustrating  time-tagged  round  trip  time  of  flight  of  the  optical  pulses  used  to 
determine  range  to  spacecraft.  The  intensity  distribution  on  the  ground  is  the  far  field  diffraction  pattern  generated 
by  the  size  and  shape  of  the  target’s  aperture.  The  satellite’s  apparent  orbital  velocity  determines  the  offset  between 
the  central  maximum  of  the  diffraction  pattern  and  the  ground  station. 

results  that  validate  the  numerical  models  for  optical  properties.  Section  7  summarizes  the  space 
qualification  testing  of  the  retroreflector  array;  Section  8  is  the  conclusion.  Three  appendices  follow 
that  cover  the  analytic  details  of  the  numerical  methods,  the  full  space  qualification  report,  and 
source  code  listings.  An  acronym  list  is  provide  at  the  end  of  the  main  text. 


2  SATELLITE  LASER  RANGING  SYSTEM  CHARACTERISTICS 

In  the  context  of  this  report,  satellite  laser  ranging  is  direct-detection  radar  in  the  optical 
wavelength  regime.  When  an  orbit  is  properly  sampled  with  a  well-calibrated  SLR  data  acquisition 
system,  residuals  and  accuracies  on  the  order  of  centimeters  are  routinely  obtained  [2]. 

Figure  1  illustrates  basic  aspects  of  the  technique.  Time-tagged  time-of-flight  differences  are 
recorded,  corrected  for  system  delays,  and  converted  to  ranges.  These  ranges  then  provide  input 
to  an  orbit  determination  model  that  is  used  to  generate  a  three-dimensional  estimate  of  the 
spacecraft’s  orbit  and  position. 
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2.1  Return  Pulse  Detection  and  System  Trade-offs 


The  number  of  photoelectrons  generated  by  an  SLR  system  is  given  by  the  laser  radar  link 
equation 

A'pe  =  VdEq  tjtGtc lrcs  (4^2)  ArtjrTZTc  .  (1) 

The  factors  in  Eq.  (1)  are: 


V D 

detector  quantum  efficiency 

O'LRCS 

laser  radar  cross  section 

E0 

transmit  energy 

R 

slant  range 

A 

wavelength 

Ar 

receiver  telescope  area 

h 

Planck’s  constant 

m 

receiver  efficiency 

c 

speed  of  light 

Ta 

one-way  atmospheric  transmission 

nr 

transmission  efficiency 

Tc 

one-way  cirrus  cloud  transmission 

Gr 

transmitter  gain 

A  requirement  of  Npe  >  10  is  a  conservative  standard  for  link  closure  and  is  used  in  this  report. 
Although  transfer  efficiencies  along  telescope  optical  paths  have  been  relatively  stable  over  the  last 
few  decades,  detector  efficiency  continues  to  improve.  Consequently,  it  is  useful  to  have  the  return 
pulse  photons  reaching  the  detector  N-,  as  a  system  figure  of  merit: 


N-, 


Npe 

VD 

Eq  (-)  m-Grcmcs  (4^2)  ArVrEIT^ 


(2) 


The  link  closure  requirement  for  photons  is  1V7  >  100.  Equation  (2)  depends  on  the  telescope 
aperture  Ar.  If  the  photon  flux  itself  is  needed  for  ground  station  trade  studies,  it  is  given  by 


flUXy  = 


Npe 

Artirtid 


(3) 


Based  on  the  ground  station  specification  in  Section  2.2,  the  results  of  this  report  are  given  in  terms 
of  Npe  and  N~t  for  a  28-cm  aperture  and  a  1-m  aperture,  as  well  as  a  13-cm  aperture  for  an  extreme 
limiting  case. 

The  factors  in  Eq.  (1)  determine  the  SLR  trade-space  and  can  be  grouped  into  four  categories: 
transfer  efficiencies,  transmitted  pulse  magnitude,  environmental/orbit  parameters,  and  geometric 
factors.  The  three  transfer  efficiencies,  tjd,  t]t-  and  tjr,  are  fixed  by  the  technology  of  the  ground 
station.  The  initial  number  of  photons,  Eq  (]£:)>  is  determined  by  the  laser  source.  The  three 
environmental/orbital  parameters  in  Eq.  (1),  i?,  Ta,  and  Tc,  are  usually  fixed  boundary  conditions 
beyond  the  experimenters’  control.  Two  of  the  geometric  factors,  Gj  and  Ar.  are  associated  with 
transmitting  and  receiving  the  laser  ranging  pulse  at  the  ground  station.  For  a  Gaussian  beam 
profile,  the  transmitter  gain  is 

Gt  =  e~2(W/*r)2  }  (4) 

where  0t  is  the  divergence  half-angle  and  0point  is  the  pointing  uncertainty.  The  remaining  geometric 
factor,  ctlrcSi  is  determined  by  the  SLR  target  and  is  a  primary  focus  of  this  report’s  analyses.  It 
is  this  factor  that  can  most  effectively  compensate  for  the  range  loss  that  goes  as  R~4.  Degnan’s 
review  [2]  provides  an  extended  treatment  of  each  of  the  parameters  in  Eq.  (1). 
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Table  1  -  Transmitter-Receiver  Configurations 


NASA/TLRS@MRC 

NRL@MRC 

FTLRS@Capraia 

Wavelength  (nm) 

532 

532  or  1064 

532 

Average  Energy  (mJ) 

100 

300 

100 

Pulse  Width  (ps) 

100  at  5  Hz 

250  at  10  Hz 

100  at  5  Hz 

Transmission  Path 

66% 

55% 

50% 

Efficiency 

Receiver  Path  Efficiency 

54% 

26%  day 
50%  night 

— 

Detector  Efficiency 

16% 

16% 

— 

Combined  Reciver  Path/ 

— 

— 

20% 

Detector  Efficency 

Full  Angle  Beam 

100 

100 

100 

Divergence  (//rad) 

Pointing  Accuracy  (arcsec) 

5.4 

1  — >•  2 

5.4 

2.2  Ground  Station  Specifications  and  LEO  Orbit  Characteristics 

To  maximize  opportunities  for  ranging  and  orbital  sampling  over  the  site,  the  NRL  LEO  retrore- 
flector  array  was  designed  to  provide  a  robust  return  margin  for  elevation  angles  down  to  the  Federal 
Aviation  Administration  (FAA)  operation  limit.  Specifically,  the  NRL  LEO  retroreflector  array  de¬ 
sign  provides  a  <tlrcs  >  104  m2  throughout  the  arc  above  an  elevation  angle  of  20°  during  a  pass, 
for  the  limiting  case  of  a  TLRS  link. 

The  link  analyses  presented  in  this  report  are  grouped  into  two  categories.  The  primary  results 
for  the  NRL  LEO  retroreflector  array  were  computed  using  a  ground  station  at  MRC  and  a  typical 
LEO  orbit,  nominally  circular,  with  an  altitude  of  1,100  km  and  a  63.5°  inclination.  The  MRC  site 
location  is  latitude  38°29'52.0770"  N;  longitude  75°22T3.4726"  W  at  30.3  m  above  sea  level.  We 
have  also  compared  link  margins  for  a  single  retroreflector  and  the  retroreflector  array  using  the 
field  transportable  laser  radar  station  (FTLRS)  and  a  TOPEX/Poseidon  orbit.  The  ground  station 
assumed  for  this  second  set  of  link  analyses  is  located  on  Capraia  Island,  latitude  43°0'  N,  longitude 
9°45/  E.  This  telescope  has  a  13-cm  aperture  and  represents  one  of  the  extreme  limiting  cases  in  the 
international  SLR  network.  TOPEX/Poseidon  satellite  and  orbit  is  of  interest  because  the  satellite 
flies  in  a  circular  orbit  at  an  altitude  of  1,330  km  with  a  66.0°  inclination.  The  spacecraft  itself 
has  a  GPS  receiver  onboard  as  well  as  a  ring  of  retroreflectors.  These  payloads  enable  independent 
orbit  determination  using  differential  GPS  as  well  as  the  SLR  data  type  [3] . 

The  MRC  link  analyses  are  based  on  two  different  transmitter-receiver  configurations  (Table  1). 
The  first  set  used  the  parameters  characteristic  of  NASA’s  TLRS.  The  second  link  was  computed 
based  on  the  NRL  Mobile  Optical  Data  Collection  Site,  which  is  anticipated  to  see  “first  light”  in 
mid-2000.  The  performance  figures  for  this  link  should  match  the  NRL  SLR  system  deployed  at 
the  Starfire  Optical  Range  (SOR)  scaled  from  SOR’s  3.5  m  to  a  1-m  aperture  telescope  [4].  FTLRS 
link  analyses  parameters  are  also  shown  in  Table  1. 

The  atmospheric  transmission  Ta  was  determined  as  a  function  of  elevation  using  MODTRAN. 
Figure  2  compares  atmospheric  transmission  as  a  function  of  elevation  at  a  midlatitude  low-elevation 
coastal  site  for  both  532  nm  and  1064  nm  wavelengths.  As  can  be  seen  from  this  graph,  transmission 
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Fig.  2  -  MODTRAN  estimates  for  atmospheric  transmission  Ta  as  a  function  of  elevation  angle  9eie v  for  a  clear  sky 
at  532  nm  (green)  and  1064  nm  (near  infrared)  for  a  typical  midlatitude  coastal  site 


at  1064  nm  is  better  than  at  532  nm.  Although  the  NASA  SLR  system  does  not  use  1064  nm,  this 
wavelength  could  be  available  in  the  NRL@MRC  1-m  system.  Furthermore,  1064  nm  is  not  visible 
to  the  naked  eye;  it  is  not  likely  to  disturb  neighboring  households  and  is  potentially  more  covert. 
However,  internal  gains  for  detectors  in  the  near-infrared  regime  are  significantly  lower  than  those 
in  the  visible  (about  a  factor  of  103).  Therefore,  signal  amplification  in  the  voltage  regime  may  be 
required  and  the  related  signal-to-noise  issues  considered.  For  all  link  analyses,  we  assume  that  no 
cirrus  clouds  are  present  and,  hence,  Tc  =  1.  All  links  can  detect  returns  in  the  daytime,  at  the 
day/night  terminator  boundary,  and  at  night. 


2.3  Target  Diffraction  Effects  and  ctlrcs 

The  typical  target  for  an  SLR  application  uses  one  or  more  retroreflectors.  A  corner  cube 
retroreflector  is  made  of  three  mutually  perpendicular  reflective  flats.  The  configuration  returns 
light  along  the  path  of  incidence,  and  therefore,  avoids  a  requirement  for  precise  orientation  of 
the  target  with  respect  to  the  ground  site.  However,  because  the  retroreflector’s  aperture  is  of 
finite  size,  the  reflected  light  is  spread  by  diffraction.  The  transmitted  pulse  travels  multiple  kilo¬ 
meters,  and  over  the  retroreflector’s  aperture  the  wave  fronts  are  indistinguishable  from  a  plane 
wave.  When  combined  with  the  return  path,  the  geometry  satisfies  the  Fraunhofer  limit.  Conse¬ 
quently,  the  spatial  distribution  of  the  return  pulse  is  the  aperture’s  far  field  diffraction  pattern 
(FFDP).  Furthermore,  the  FFDP  follows  standard  diffraction  scaling  in  which  a  smaller  aperture 
retroreflector  generates  a  broader  FFDP.  The  FFDP’s  breadth  is  of  particular  importance  since,  in 
combination  with  the  satellite  velocity  aberration,  it  plays  a  significant  role  in  overall  SLR  system 
performance. 
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Fig.  3  -  Relation  between  k  and  the  angles  9X  and  9y. 


In  Eq.  (1),  the  SLR  target’s  diffraction  characteristics  are  carried  by  ctlrcs •  Denoting  the  FFDP 
by  a(kx,ky),  the  relationship  between  the  two  quantities  [5,  6]  is 

47T 

0LRCs{kx,  ky)  =  -^\a(kx,  fc^)!2  ,  (5) 

where 

7  27T 

kx  =  —  sin  9X  (6) 

and 

7  2tt  .  A 

ky  =  Tsm^  '  (7) 

The  angles  6X  and  9y  define  the  observation  direction.  If  an  FFDP  is  observed  at  x0  and  y0  from  a 
range  R,  then  sin  9X  =  x0/R  and  sin  9y  =  ya/R,  respectively.  One  is  equally  free  to  use  a  position 
and  a  distance,  but  for  general  work  it  is  more  convenient  to  use  angles.  The  results  of  this  report 
are  given  in  terms  of  angular  measure.  Figure  3  shows  the  relation  between  k  and  the  angles  6X  and 
9y.  Equation  (5)  assumes  perfect  reflectivity.  A  more  realistic  description  of  a  real  retroreflector  is 

^lrcs(^,  ky)  =  \a{kx,  ky)\2  ,  (8) 

where  p  is  a  reflectivity.  In  the  computations  supporting  this  report,  p  =  0.75,  which  is  a  conser¬ 
vative  estimate  that  includes  manufacturing  errors  as  well  the  reflectivity  loss  itself. 


2.4  Satellite  Velocity  Aberration 

The  terms  “velocity  aberration,”  “relativistic  effect,”  and  “Bradley  effect”  all  refer  to  the  same 
physical  process.  When  light  signals  are  observed  from  a  system  moving  with  respect  to  the  source, 
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Fig.  4  -  The  NRL  LEO  retroreflector  array  using  22  individual  retro  reflectors  on  the  surface  of  hemisphere.  This 
array  has  a  FOV  of  108°  and  orLRcs  >  104  m2  for  Qeiev  >  20°. 

their  angular  relationships  change.  In  the  case  of  SLR,  the  satellite’s  motion  relative  to  the  ground 
station  causes  the  part  of  the  intensity  distribution  from  the  return  pulse  that  is  observed  at  the 
ground  station  to  be  offset  with  respect  to  the  central  maximum  of  the  diffraction  pattern.  In  SLR 
parlance,  the  most  common  terms  for  this  offset  are  satellite  velocity  aberration  or  the  Bradley 
effect  [5].  Figure  1  shows  schematically  the  combined  effects  of  the  diffracted  return  pulse  and 
the  velocity  aberration.  The  velocity  aberration  expressed  as  a  displacement  vector  in  the  far  field 
diffraction  pattern  (or  ctlrcs)  *s  giyen  by 

£=(f)(^)  ■  <9> 

where  vj_  is  the  satellite’s  apparent  perpendicular  velocity.  The  displacement  vector  k  selects  which 
value  of  the  instantaneous  olrcs  from  Eq.  (8)  is  used  in  Eq.  (1)  at  each  point  in  time  during  a 
satellite’s  pass  over  a  ground  station. 

3  NRL  LEO  RETROREFLECTOR  ARRAY 

The  NRL  LEO  array  comprises  22  1-cm  retroreflectors  specifically  oriented  in  a  82-mm  diameter 
hemisphere  held  in  an  aluminum  mount  43  mm  in  height  yielding  a  total  mass,  with  bolts,  of 
221  grams.  The  flight-qualified  retroreflectors  have  0.004-in.  bevels.  The  retroreflectors  are  inset 
slightly,  which  results  in  some  vignetting  (i.e.,  c  ^  0  in  Eq.  (A8)).  The  qualification  unit  has  0.008- 
in.  bevel  retroreflectors  and  is  the  basis  for  the  analyses  and  experiments  in  this  report.  Figure  4 
is  a  photograph  of  the  array. 

The  placement  of  the  retroreflectors  in  an  array  is  determined  primarily  by  its  effect  on  the 
spatial  distribution  of  return  pulse  optical  energy,  rather  than  strictly  mechanical  properties.  In 
other  words,  the  location  of  the  retroreflectors  should  be  likened  to  placement  of  the  elements  of 
an  RF  antenna.  If  the  retroreflectors  are  placed  incorrectly,  it  is  possible  for  the  array’s  diffraction 
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pattern  to  combine  with  the  velocity  aberration  and  generate  a  null  at  the  ground  station.  The 
NRL  LEO  array’s  multiple  retroreflectors  are  specifically  placed  to  provide  a  field  of  view  (FOV) 
of  108°,  while  neither  generating  such  a  null  nor  impacting  timing  (see  Section  5.3).  Figure  5  is  a 
mechanical  drawing  of  the  array.  The  retroreflectors  are  in  three  rings.  The  first  ring  is  at  16°  with 
four  retroreflectors  placed  at  0°,  90°,  180°,  and  270°  in  azimuth.  The  second  is  at  32°  with  eight 
retroreflectors  placed  at  22.5°,  67.5°,  112.5°,  157.5°,  202.5°,  247.5°,  292.5°,  and  337.5°  in  azimuth. 
The  third  is  at  48°  with  10  retroreflectors  placed  at  18°,  54°,  90°,  126°,  162°,  198°,  234°,  270°, 
306°,  and  342°  in  azimuth. 

To  compute  predicted  on-orbit  SLR  performance  for  this  array,  both  far  field  diffraction  ef¬ 
fects  and  velocity  aberration  must  be  taken  into  account.  Our  numerical  approach  required,  as  a 
building  block,  that  we  first  develop  a  method  to  predict  the  on-orbit  SLR  performance  of  a  single 
retroreflector.  Consequently,  we  can  compare  the  predicted  performance  between  the  NRL  LEO 
array  with  that  of  a  single  retroreflector.  To  foreshadow  Section  5.2,  the  returns  from  the  array’s 
multiple  retroreflectors  for  a  given  pulse  mitigates  nulls  of  a  single  retroreflector  ctlrcs  and  permit 
a  much  larger  FOV. 

4  NUMERICAL  COMPUTATION  OF  SATELLITE  LASER  RANGING 
PERFORMANCE 

We  have  developed  a  set  of  numerical  techniques  to  predict  on-orbit  SLR  performance.  The 
three  basic  problems  addressed  are: 

1.  computing  the  ctlrcs  of  a  single  retroreflector; 

2.  combining  several  single  retroreflector  ctlrcs  values  to  determine  an  array’s  overall  <7LRCS; 
and 

3.  using  orbit  dynamic  data  with  the  <tlrcs  values  to  compute  the  actual  link  analysis. 

This  section  sketches  the  physical  basis  for  the  numerical  methods  and  the  relative  impact  of 
different  SLR  properties  on  performance.  The  underlying  analytic  expressions  and  numerical  im¬ 
plementation  are  presented  in  detail  in  Appendix  A  and  Appendix  B  respectively.  The  experimental 
validation  of  results  from  these  routines  is  discussed  in  Section  6. 

4.1  Single  Circular  Retroreflector  Far  Field  Diffraction  Patterns  and  ctlrcs 

This  section  discusses  spatial  distribution  of  a  pulse  return  from  a  single  retroreflector  caused  by 
diffraction.  The  section  covers  the  results  from  our  numerical  computation  of  ctlrcs  for  four  cases 
beginning  with  the  highest  symmetry:  normal  incidence  without  bevels;  tilted  incidence  without 
bevels;  normal  incidence  with  bevels;  and  tilted  incidence  with  bevels.  Each  case  serves  as  a  limiting 
test  for  the  numerical  routines  of  the  succeeding  cases. 


4-1-1  Normal  Incidence 

When  a  plane  wave  pulse  hits  the  retroreflector,  it  is  exactly  reversed  in  direction  along  the 
path  of  incidence.  As  discussed  in  Section  2.3,  the  retroreflector’s  finite  aperture  modifies  the 
incidence  plane  wave  pulse  and  returns  a  diffracted  pulse.  At  the  ground  station,  the  return 
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pulse  has  the  spatial  distribution  of  the  aperture’s  far  field  diffraction  pattern.  A  single  circular 
retroreflector  without  bevel  losses  observed  at  normal  incidence  provides  a  geometry  with  sufficiently 
high  symmetry  that  an  analytic  expression  exists  for  the  FFDP.  The  FFDP  is  the  Airy  function 
[7,  8];  hence,  from  Eq.  (5), 


^LRCsW  = 


4  ttA2 


(  2Ji(rfc)\  2 
V  rk  ) 


(10) 


where  A  is  the  area  of  the  retroreflector,  r  is  the  retroreflector  radius,  and 


k  =  sin  9 

A 


At  normal  incidence,  ctlrcs  has  azimuthal  symmetry  and  Eq.  (10)  is  a  function  of  the  single 
magnitude  variable  k  only.  Note  that  in  Eq.  (11),  6  is  the  interior  angle  between  the  +Z-axis  and 
the  observation  direction  in  Fig.  3.  It  is  also  more  convenient  to  plot  ctlrcs  in  terms  the  equivalent 
0 .  Even  for  a  ctlrcs  f°rm  that  lacks  azimuthal  symmetry,  kx  and  ky  arguments  can  be  specified  in 
terms  of  a  pair  of  equivalent  angles  6X  and  0y. 

Figure  6  shows  ctlrcs  as  function  of  k  (a  radial  slice  in  kx ,  ky)  for  2.54-cm  and  1-cm  aperture 
sizes.  Both  curves  were  computed  numerically  with  the  FFDP  routines  described  in  Appendix  A.l 
and  agree  with  Eq.  (10)  to  machine  precision.  The  correct  geometry  required  for  normal  incidence 
can  occur  only  when  the  satellite  is  directly  above  the  ground  station  and  the  retroreflector  normal 
is  aligned  with  the  nadir  direction.  The  velocity  aberration  magnitude  varies  approximately  over 
the  range  39  to  49  /irad  for  a  1,100-km  circular  orbit  (typical  LEO  spacecraft)  and  over  38  to 
48  //rad  for  a  1,330-km  circular  orbit  (TOPEX/Poseidon).  Thus,  as  can  be  seen  in  Fig.  6,  although 
ctlrcs  for  the  2.54-cm  retroreflector  is  significantly  larger  at  k  =  0,  the  broader  FFDP  of  the  1-cm 
retroreflector  can  have  a  greater  LRCS  over  the  critical  regions  and  certainly  has  greater  constancy 
(less  variation)  throughout  the  pattern.  Figure  7  is  a  contour  plot  of  ctlrcs  for  a  1-cm  circular 
retroreflector  as  a  function  of  kx  and  ky  in  terms  of  //rad. 

While  Fig.  7  shows  ctlrcs  over  a  sizable  angular  region  of  the  far  field,  the  observation  of  return 
pulse  is  made  at  a  single  kXl  ky  point.  The  velocity  aberration  given  by  Eq.  (9)  selects  this  single  kx , 
ky  point  and,  consequently,  the  ctlrcs  value  used  in  Eq.  (1)  during  the  link  margin  computation. 
To  provide  a  scale  length,  for  an  1,100-km  orbit,  the  corner  points  of  Fig.  7  at  ±240  ^radians 
correspond  to  ±0.264-km  distances. 


4-1.2  Off-axis  Incidence 

At  all  other  points  in  an  orbital  pass,  the  retroreflector  tips  away  from  normal  incidence,  so  its 
projected  aperture  changes  in  shape  and  decreases  in  area.  Figure  8  shows  how  both  the  shape 
of  the  projected  aperture  and  the  area  changes  as  a  function  of  incidence  angle  &inci.  Note  that 
the  projected  aperture  decreases  in  size  along  both  axes.  The  decrease  in  projected  aperture  size 
increases  the  angular  extent  of  the  FFDP  proportionally  and  decreases  its  overall  magnitude.  The 
change  in  aperture  shape  eliminates  azimuthal  symmetry  so  an  analytic  closed-form  expression 
for  the  FFDP  (e.g.,  Airy  function)  is  no  longer  possible.  Appendix  A.l  describes  a  method  to 
numerically  compute  the  FFDP. 

Figure  9  shows  ctlrcs  of  a  tilted  1-cm  circular  retroreflector  computed  numerically  without 
bevel  losses.  The  azimuthal  symmetry  present  for  normal  incidence  has  been  reduced  to  a  two-fold 
rotational  symmetry  as  seen  from  the  ellipticity  of  the  contours.  The  central  maximum  around 
kx  —  ky  =  0  remains  smoothly  rounded.  For  observations  made  away  from  the  central  maximum, 
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Fig.  8  -  (Top)  Retroreflector  aperture  shape:  (left  to  right)  Normal  incidence;  20°  from  normal  incidence;  40°  from 
normal  incidence.  (Bottom)  Normalized  retroreflector  aperture  area  as  a  function  of  incidence  angle  0inci .  FOV 
cutoff  is  Qind  =  58.579°,  which  places  a  significant  limitation  on  orbit  sampling.  See  Section  5.2. 


however,  there  is  a  sizable  variation  in  ctlrcs  as  one  holds  the  magnitude  of  k  fixed  and  sweeps 
around  in  kx,  ky.  Physically,  this  corresponds  to  tilting  the  retroreflector  axis  to  a  fixed  orienta¬ 
tion  while  sweeping  the  observation  direction  around  the  surface  of  a  cone.  Figure  10  shows  this 
variation  in  ctlrcs  f°r  a  10°  tilt  from  normal  incidence  and  |k|  equivalent  to  50  yurad  as  a  function 
of  observation  direction  angle.  Note  that  50  /irad  is  approximately  the  maximum  LEO  velocity 
aberration.  For  smaller  velocity  aberrations,  these  variations  in  ulrcs  are  correspondingly  smaller 
as  well. 


4-1-3  Bevel  Losses 

Beveled  edges  are  required  at  the  reflecting  face  plane  intersections  during  the  polishing  of 
real  retroreflector  cubes.  Adding  bevels  introduces  narrow  loss  regions  within  the  aperture.  More 
general  versions  of  our  numerical  simulations  include  these  losses  for  both  normal  incidence  and 
off-normal  incidence  cases  (Sections  A.  1.3  and  A.  1.4).  At  normal  incidence,  the  three  bevel  loss 
regions  and  their  reflections  generate  a  six-fold  rotation  symmetry  in  diffraction  pattern  and  clrcs- 
Since  the  bevels  are  relatively  narrow,  the  six-fold  symmetry  is  more  pronounced  at  larger  values  of 
|k|  beyond  the  first  diffraction  minimum.  Figure  11  shows  clrcs  normal  incidence  of  a  circular 
retroreflector  computed  numerically  with  losses  from  0.008-in.  bevels.  Figure  12  show  the  variation 
in  0LRCS  as  one  sweeps  around  in  kX)  ky  for  0.008-in.  bevels  at  |ic|  equivalent  to  50  /xrad.  The 
variation  in  Fig.  12  is  only  0.25%. 

Comparing  Figs.  10  and  12  shows  that  the  effect  of  tilting  the  retroreflector  is  much  more 
significant  than  adding  the  bevels.  Although  our  numerical  computation  treats  the  bevel’s  effect 
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Fig.  9  -  (Tlrcs  of  a  1-cm  circular  retroreflector  as  a  function  of  kx  and  ky  given  in  terms  of  /irad,  with  10°  tilt. 
Ellipticity  of  the  contours  is  the  result. 


Fig.  10  -  Variation  in  o'lrcs  at  10°  off-axis  from  normal  incidence  of  a  1-cm  circular  retroreflector  at  |k|  equivalent  to 
50  /irad.  The  variation  in  oxrcs  from  its  average  caused  by  a  fixed  10°  tilt  by  a  rotation  of  the  observation  direction 
is  42%.  At  the  instant  of  pulse  reflection,  a  satellite’s  apparent  perpendicular  velocity  determines  the  observation 
direction. 
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equivalent  to  50  prad  with 


exactly,  to  first  order,  the  bevels  could  also  be  considered  as  simply  a  scalar  loss  of  reflecting  area 
combined  with  a  slightly  reduced  FOV.  The  combination  of  both  retroreflector  tilt  and  bevel  losses 
requires  a  slightly  more  complex  numerical  simulation,  but  it  does  not  introduce  any  qualitatively 
new  features  to  the  FFDP  or  <tlrcs- 


4.2  Retroreflector  Array  Far  Field  Diffraction  Patterns  and  ctlrcs 


The  return  pulse  is  the  sum  of  the  contributions  from  several  retroreflectors  in  the  array.  The 
total  FFDP  a  for  the  array  of  L  retroreflectors  at  a  given  instant  of  time  is 

a  =  f>e‘“‘  ,  (12) 


where  the  eqs  are  the  phase  angles  for  each  retroreflector  for  a  given  incidence  direction.  From 
Eq.  (8),  the  instantaneous  LRCS  is 


I--*  I 

O'LRCS  =  P-jj\aa  I  • 


Upon  substitution: 


tf'LRCS 


■  '5 1(5,  *-•“)  (S5-™) 

=  EEW"”-”"’  • 


During  a  brief  interval  r,  variations  of  phase  angle  relationships  between  individual  retroreflec¬ 
tors  can  occur.  Provided  the  statistical  correlation  between  am  and  an  over  r  is  small,  one  is  free 
to  treat  resulting  phase  factors  as  essentially  random  fluctuations.  Consequently,  a  time-averaged 
LRCS  will  simplify  to: 


^LRCS 


P~  TZ  V  Qmfl"  fT  e^a(t)m~a(t^dt 

A2  m=\  n=l  T 


4tt  i~ 

P  \2  I0*1 


(16) 

(17) 


For  SLR  measurements  with  precisions  significantly  tighter  than  ±1  cm,  r  may  be  become  small 
enough  that  the  phase  integral  cannot  be  considered  as  a  Kronecker  delta  function  Smn.  For  the 
LEO  application,  however,  using  a  random  phase  approximation  is  sufficient,  and  the  array’s  LRCS 
can  be  evaluated  using  the  incoherent  sum  of  the  contributions  in  Eq.  (17). 

Equation  (17)  is  a  convenient  form  for  numerical  evaluation.  The  routines  developed  and  tested 
in  the  single  retroreflector  can  be  used  directly  to  compute  the  am s.  The  primary  addition  is 
a  routine  that  combines  an  incidence  direction  with  respect  to  the  array’s  axis  with  individual 
retroreflector  directions  to  determine  which  retroreflectors  are  illuminated  within  their  FOV  limits 
and  at  what  angles.  It  is  then  straightforward  to  compute  the  relevant  aTOs  and  to  sum  to  obtain 
<tlrcs.  Figure  13  shows  ctlrcs  of  the  array  (using  retroreflectors  with  0.008-in.  bevels)  as  computed 
numerically  when  viewed  along  the  array’s  axis.  (Note  that  array  ctlrcs  values  used  in  Section  5 
are  time-averaged.  Unless  there  is  a  specific  reason  to  do  otherwise,  explicit  notation  of  the  time- 
averaging  is  suppressed  for  the  balance  of  this  report.) 
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Fig.  13  -  NRL  LEO  retroreflector  array  <tlrcs  with  0.008-in.  bevels  and  on-axis  incidence 


4.3  Computation  of  Orbital  Performance 

The  local  geometry  at  the  satellite  and  apparent  velocity  for  a  given  ground  station  must  be 
used  to  drive  the  ctlrcs  roufines  to  compute  the  performance  expected  for  a  either  the  NRL  LEO 
retroreflector  array  or  a  single  retroreflector  on  orbit.  The  range  vector,  apparent  velocity  vector, 
and  local  elevation  angle  from  nadir  are  sufficient  to  determine  the  kXl  ky  point  and  the  form  of 
<7lr.cs  that  will  occur  at  a  given  point  in  time  during  a  satellite’s  orbital  pass.  (See  Section  A. 2  for 
the  detailed  development.)  The  SLR  performance  can  then  determined  with  Eq.  (1)  at  each  point 
in  time  during  a  pass  for  which  the  satellite’s  attitude  and  apparent  velocity  are  known. 

Computation  of  a  spacecraft’s  orbital  dynamics  sufficient  to  determine  the  instantaneous  atti¬ 
tude  and  apparent  velocity  is  a  separate  and  specialized  task.  The  ephemeris  and  attitude  data 
used  in  this  report  were  supplied  by  NRL’s  Electro-Optics  Technology  Section  and  Mathematics  & 
Orbit  Dynamics  Section  (Codes  8123  and  8233). 

5  PREDICTED  ORBITAL  PERFORMANCE 

5.1  NRL  LEO  Array  Performance  in  Orbit 

As  previously  discussed,  the  local  geometry  and  velocity  aberration  combine  to  determine  olrcs 
at  any  given  instant.  Equation  (1)  can  then  be  used  to  determine  Npe  as  a  function  of  time  during 
a  pass  over  the  ground  site.  For  a  nadir-facing  array,  the  local  angle  from  nadir  is  the  angle  between 
the  incidence  direction  and  the  array’s  axis.  The  apparent  velocity  vector  is  similarly  common  to 
the  entire  array.  As  noted  in  Section  4.2,  this  information  is  first  converted  to  the  coordinate  system 
of  each  illuminated  retroreflector  in  the  array  and  then  used  with  Eq.  (17)  to  compute  ctlrcs. 
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Fig.  14  -  ctlrcs  of  the  22-element  NRL  LEO  retroreflector  array  at  532  and  1064  nm  (top)  and  $etev  (bottom)  as  a 
function  of  time;  this  is  a  near-zenith  pass  over  MRC  with  0eiev  =  88°  at  PC  A 


Figures  14-23  show  the  NRL  LEO  retroreflector  array’s  performance  on  orbit  from  the  MRC 
ground  site.  These  computations  are  based  on  the  link  parameters  for  the  TLRS  and  the  NRL@MRC 
system  described  in  Section  2.2,  as  well  as  the  MODTRAN  atmospheric  transmission  values  as  a 
function  of  pass  elevation  shown  in  Fig.  2.  We  analyzed  the  link  for  two  typical  passes  over  the 
MRC  site:  one  that  has  a  near-zenith  elevation  angle  of  88°  at  point-of-closest  approach  (PCA), 
and  a  second,  much  lower,  pass  with  6eiev  =  39°  at  PCA.  This  comparison  determined  the  impact 
of  low  pass  geometries  and  velocity  aberration  on  link  closure  using  the  retroreflector  array. 

Figure  14  shows  ctlrcs  at  532  and  1064  nm  vs  time  as  well  as  6eiev  vs  time  for  the  near-zenith 
pass  over  MRC.  Figure  15  shows  the  equivalent  three  curves  of  532  and  1064  nm  ctlrcs  and  Qeiev 
vs  a  common  time  axis  for  the  low  elevation  off-zenith  pass.  As  can  be  seen  from  Fig.  14,  although 
there  is  a  drop  at  PCA,  ctlrcs  remains  greater  than  104  m2  throughout  the  entire  pass.  For  the  low 
elevation  pass,  Fig.  15  shows  ctlrcs  >  104  m2  occurring  even  at  the  lowest  elevation  angles.  The 
asymmetry  apparent  in  ctlrcs  during  the  pass  is  produced  by  variations  in  satellite  attitude  and 
velocity  aberrations.  These  effects  are  more  pronounced  than  in  a  zenith  pass.  The  ctlrcs  values 
plotted  in  Figs.  14  and  15  were  used  in  Eqs.  (1)  and  (2)  for  the  link  analyses  shown  in  Figs.  16-23. 


In  Figs.  16  and  17,  the  number  of  return  photons  N 7  generated  at  532  nm  reaching  the  TLRS 
telescope’s  detector  are  shown  during  high  and  low  passes,  respectively.  Figures  18  and  19  show 
the  number  of  photoelectrons  Npe  generated  as  a  function  of  time  for  these  passes  at  532  nm  using 
a  photomultiplier  tube  with  a  quantum  efficiency  of  16%  and  an  internal  gain  of  106,  which  is 
standard  technology.  Figure  19  clearly  shows  that  the  NRL  LEO  array  can  close  the  link  with 
Npe  >  10  at  the  low  elevations,  even  with  the  28-cm  TLRS  telescope. 
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Fig.  15  -  «tlrcs  of  the  22-element  NRL  LEO  retroreflector  array  at  532  and  1064  nm  (top)  and  9eiev  (bottom) 
function  of  time;  this  is  a  low-elevation  pass  over  MRC  with  0etev  =  39°  at  PC  A 
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Fig.  16  -  N-y  at  532  nm  for  a  28-cm  diameter  aperture  TLRS  telescope  and  the  NRL  LEO  retroreflector  array  at 
MRC  (top)  and  0ejet?  (bottom)  as  a  function  of  time  during  a  near-zenith  pass 
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Fig.  17  -  N*y  at  532  nm  for  a  28-cm  diameter  aperture  TLRS  telescope  and  the  NRL  LEO  retrorefiector  array  at 
MRC  (top)  and  0eiev  (bottom)  as  a  function  of  time  during  a  low  pass 
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Fig.  18  -  Npe  at  532  nm  for  a  28-cm  diameter  aperture  TLRS  telescope  and  the  NRL  LEO  retrorefiector  array  at 
MRC  (top)  and  9eiev  (bottom)  as  a  function  of  time  during  a  near-zenith  pass 
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Fig.  19  -  Npe  at  532  nm  for  a  28-cm  diameter  aperture  TLRS  telescope  and  the  NRL  LEO  retroreflector  array  at 
MRC  (top)  and  6euv  (bottom)  as  a  function  of  time  during  a  low  pass 


For  the  NRL  mobile  meter  class  system  to  be  situated  at  MRC,  the  link  closures  for  day  and 
night  at  532  nm  are  shown  in  Figs.  20  and  21,  again  for  high  and  low  passes.  Because  this  telescope 
system  will  be  under  NRL  control,  the  daylight  filter  can  be  added  or  removed  as  required.  Hence, 
the  number  of  received  photoelectrons  can  be  increased  during  terminator  and  at  night.  The 
advantage  of  this  clearly  more  robust  link  is  that  more  unambiguous  yield  is  possible  throughout  a 
pass,  particularly  at  the  lower  elevation  angles  and  daytime  operation.  This  yield  is  sufficiently  high 
that  a  postdetector  amplifier  is  not  required,  and  the  associated  amplifier  noise  can  be  eliminated. 
Thus,  orbital  sampling  over  the  site  is  maximized  and  precision  is  improved. 

Transmission  through  the  atmosphere  at  1064  nm  is  better  than  at  532  nm  as  previously  dis¬ 
cussed  (see  Fig.  2).  Since  near-infrared  detector  technology  can  be  expected  to  show  significant 
changes,  N 7  is  a  more  appropriate  figure  of  merit.  Figures  22  and  23  show  the  N 1  curves  for  day¬ 
time  1064  nm  observations  at  MRC.  As  discussed,  comparable  internal  gains  are  not  available  for 
detectors  in  the  near-infrared.  Furthermore,  signal-to-noise  issues  become  significant  at  1064  nm. 
However,  quantum  efficiencies  can  be  higher  (20%  to  40%).  Operation  at  1064  nm  does  have  the 
additional  inherent  advantage  of  being  well  outside  the  range  of  human  vision. 

5.2  NRL  LEO  Array  and  Single  Retroreflector  Performance  Comparison 

To  directly  compare  of  performance  between  a  single  retroreflector  and  the  NRL  LEO  array,  we 
considered  a  specific  pass  with  TOPEX/Poseidon  orbital  parameters  as  the  test  case.  Note  that 
a  different  midlatitude  ground  station  is  assumed.  The  links  are  predicted  for  a  FTLRS  (13-cm 
telescope  aperture)  located  on  Capraia.  The  major  result  from  the  standpoint  of  orbit  sampling 
capability  is  apparent  directly  from  the  <jlrcs  plots  shown  in  Figs.  24  and  25.  The  NRL  LEO 
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Fig.  20  -  Npc  at  532  nm  for  a  100-cm  diameter  aperture  telescope  and  the  NRL  LEO  retroreflector  array  at  MRC 
(top)  and  deiev  (bottom)  as  a  function  of  time  during  a  near-zenith  pass 
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Fig.  21  -  NPe  at  532  nm  for  a  100-cm  diameter  aperture  telescope  and  the  NRL  LEO  retroreflector  array  at  MRC 
(top)  and  delev  (bottom)  as  a  function  of  time  during  a  low  pass 
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Fig.  22  -  N-i  at  532  nm  and  1064  nm  for  a  100-cm  diameter  aperture  telescope  and  the  NRL  LEO  retroreflector  array 
at  MRC  (top)  and  0eiev  (bottom)  as  a  function  of  time  during  a  near-zenith  pass 


Fig.  23  -  N 7  at  532  nm  and  1064  nm  for  a  100-cm  diameter  aperture  telescope  and  the  NRL  LEO  retroreflector  array 
at  MRC  (top)  and  0eiev  (bottom)  as  a  function  of  time  during  a  low  pass 
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Fig.  24  -  <tlrcs  at  532  nm  for  a  1-cm  circular  retroreflector  (top)  and  9eiev  (bottom)  as  a  function  of  time.  For 
Pass  44  at  PCA,  96iev  =  89.38°  so  this  is  effectively  a  zenith  pass.  Note  the  impact  of  retroreflector  FOV  on  <tlrcs 


array  prevents  the  drop-out  in  orbit  sampling  that  occurs  with  the  single  retroreflector  whenever 
its  narrow  FOV  drives  o*lrcs  to  zero.  For  completeness,  iVpe  as  function  of  time  for  a  single  1- 
cm  retroreflector  with  0.008-in.  bevels  and  the  NRL  LEO  array  are  shown  in  Figs.  26  and  27. 
The  performance  is  evaluated  based  on  5.4  arcsec  pointing  uncertainty  and  transmitter/receiver 
parameters  described  in  Gilbreath  et  al.  [9].  These  results  clearly  show  how  the  array  closes  the 
link  throughout  a  pass  by  opening  up  the  FOV  well  beyond  the  57°  window  that  constrains  a  single 
retroreflector. 


5.3  Optical  Phase  Center  and  Timing  Precision 

The  orbital  estimation  precision  using  the  SLR  data  type  is  on  the  order  of  a  ±1  cm  rms. 
Therefore,  knowledge  of  the  offset  of  the  retroreflector’s  location  with  respect  to  the  spacecraft 
center  of  mass  (COM)  is  important.  This  fixed  offset  is  called  the  range  correction  and  exists 
whether  a  single  retroreflector  or  an  array  is  used. 

In  the  case  of  an  array,  the  timing  precision  of  the  return  signal  can  also  be  affected  if  more 
than  one  retroreflector  is  illuminated  by  a  pulse.  The  placement  of  retroreflectors  defines  an  optical 
phase  center  (OPC)  as  the  common  origin  of  all  the  retroreflector  outward  normals.  To  the  extent 
the  array  retroreflectors  lack  a  common  OPC,  a  direction-dependent  offset  will  be  introduced  into 
the  system.  Consequently,  to  achieve  centimeter-level  precision  and  accuracy,  the  retroreflectors  in 
a  multi-element  array  must  be  aligned  with  respect  to  a  comparable  tolerance  to  a  common  surface. 

The  OPC  of  the  NRL  LEO  array  is  centered  on  the  vertical  axis  6  mm  inside  the  array.  Cen¬ 
tered  about  this  phase  center  is  a  surface  of  reflection.  The  advantage  of  a  hemisphere  of  many 
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Fig.  27  -  Npe  at  532  nm  for  a  13-cm  diameter  aperture  FTLRS  telescope  and  the  NRL  LEO  retroreflector  array  at 
Capraia  (top)  and  Oeiev  (bottom)  as  a  function  of  time  during  Pass  44;  a  5.4  arcsec  pointing  precision  is  assumed 


small  reflectors  is  that  for  any  incidence  direction  there  is  a  retroreflector  in  close  proximity.  This 
condition  allows  a  relatively  smooth  surface  of  reflection  overall.  The  radius  of  this  surface  was 
determined  by  examining  585  angles  and  calculating  the  point  on  this  surface  for  each  angle.  The 
surface  can  be  approximated  by  a  sphere  with  a  radius  of  2.505  cm.  The  standard  deviation  from 
this  sphere  is  only  0.086  cm.  This  radius  is  the  largest  contributor  to  the  range  correction. 

The  range  correction  also  must  consider  any  pulse  broadening  due  to  the  multiple  retroreflectors. 
Using  small  retroreflectors  allows  for  a  very  compact  arrangement.  The  small  size  limits  the  possible 
pulse  broadening.  The  broadening  results  from  retroreflectors  at  different  locations  on  the  array 
having  a  delayed  return.  In  the  nadir  direction,  for  example,  the  ring  of  4  retroreflectors  in  the 
center  is  closer  to  the  ground  than  is  the  middle  ring  of  8  retroreflectors.  The  time  delay  is  66  ps. 
The  return  from  the  middle  ring  is  smaller  because  of  the  larger  angle.  Combining  the  returns 
from  the  two  rings  broadens  the  return  pulse.  The  ring  of  4  has  an  average  <tlrcs  °f  52, 000  m2 
and  the  ring  of  8  has  an  average  cjlrcs  of  16,000  m2,  for  nadir.  Figure  28  shows  the  incoherent 
sum  of  the  two  sets  of  returns  for  an  incident  pulse  with  a  full  width  at  half  maximum  (FWHM) 
of  100  ps  (typical  for  TLRS).  The  return  pulse  is  broadened  to  117  ps  FWHM  with  a  7  ps  (1  mm) 
delay  in  the  peak  amplitude.  The  extra  17  ps  of  width  contributes  to  the  range  error.  A  1-mm 
two  way  range  correction  is  required  to  correct  for  pulse  broadening.  Longer  pulses  typical  of  the 
NRL@MRC  system  (250  ps)  would  have  smaller  effects. 

The  range  correction  error  for  the  time  of  flight  is  the  combination  of  the  size  of  the  sphere 
error  and  the  pulse  broadening.  The  increase  in  the  total  standard  deviation  error  of  the  two  is 
±0.86  mm  and  ±2.6  mm.  Using  Gaussian  statistics  results  in  a  total  error  of  ±2.74  mm.  This 
analysis  shows  that  the  timing  error  of  the  array  is  within  the  ±1  cm  specification.  A  more  accurate 
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Fig.  28  -  Array  pulse  timing  effects 


error  analysis  for  millimeter-level  range  precision  will  require  detailed  information,  including  the 
nonlinear  behavior  of  the  receiver  detector  and  electronics. 

It  should  be  noted  that  multiple  retroreflectors  interfere  due  to  the  coherence  of  the  laser  beam. 
This  effect  is  termed  satellite  scintillation.  An  incoherent  sum  over  the  array’s  retroreflectors 
generates  the  average  return  signal,  but  a  coherent  sum  is  required  for  determining  the  standard 
deviation  of  the  return  signal.  Having  more  than  two  retroreflectors  in  the  FOV  reduces  the 
deviation.  Since  the  location  of  each  retroreflector  at  a  given  moment  in  time  is  not  known, 
random  phases  are  assigned  to  each  retroreflector.  This  is  equivalent  to  moving  each  retroreflector 
less  than  1  fj,m.  This  was  done  analytically  for  the  clrcs  at  each  point  in  the  pass  and  was  repeated 
for  100  random  phases  for  each  retroreflector.  The  resulting  standard  deviation  of  the  ctlrcs  was 
less  than  12%  of  the  average  return  signal  for  all  855  samples  of  Pass  44  over  Capraia  (Fig.  29). 
This  deviation  is  small  compared  to  atmospheric  scintillation  [10]. 


5.4  Further  Comparison  of  Single  Retroreflectors  and  Retroreflector  Arrays 

Using  an  array  eases  the  minimum  pass  geometry  requirements  significantly.  As  has  been  shown, 
an  array  can  provide  viable  returns  at  much  lower  elevation  angles.  The  quantitative  improvement 
in  available  passes  can  be  determined  by  sorting  all  pass  geometries  from  the  10-day  repeat  cycle 
over  ground  stations  in  the  Washington,  D.C.  area  by  $eiev  at  PCA.  Figure  30  shows  the  number  of 
passes  that  exceed  a  given  value  of  0eiev  at  PCA  as  a  function  0eiev  for  the  repeat  cycle.  From  Fig.  24, 
the  single  retroreflector  clrcs  places  lower  limit  of  9eiev  >  40°,  assuming  perfect  atmospheric  seeing. 
Scintillation  and  other  effects  will  likely  restrict  sampling  to  9eiev  >  50°.  For  an  array,  the  additional 
tracking  passes  as  low  as  the  20°  FAA  limit  can  used  and,  consequently,  a  better  orbit  solution 
becomes  available. 
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Fig.  29  -  Normalized  amplitude  fluctuations  of  the  NRL  LEO  array;  rms  variation  in  average  return  signal  is  12% 


0elev  ot  PCA  (deg) 

Fig.  30  -  Number  of  passes  vs  elevation  angle  at  PCA  over  a  Washington,  D.C.  area  ground  station  for  a  typical 
1,100-km  LEO  orbit  during  a  10-day  ground  repeat  track;  dotted  lines  at  20°  and  50°  are  conservative  observation 
limits  for  the  NRL  LEO  array  and  single  1-cm  retroreflector  respectively 
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Fig.  31  -  Optical  bench  layout  for  far-field  characterization  of  the  retroreflector  array 


6  EXPERIMENTAL  DETERMINATION  OF  <jlrcs 


The  NRL  LEO  retroreflector  array’s  ctlrcs  was  measured  in  the  laboratory  to  validate  the 
numerical  computations.  This  subsection  describes  the  experimental  configuration  and  sensor  cal¬ 
ibration  procedure.  Then,  the  array’s  instantaneous  intensity  distribution  is  compared  to  the 
time-averaged  distribution.  Finally,  differences  between  experimental  and  numerical  results  are 
discussed,  with  particular  attention  to  the  increased  energy  directed  into  sidelobes  from  retroreflec¬ 
tor  imperfections. 


6.1  Procedure 

The  performance  of  the  NRL  LEO  retroreflector  array  was  experimentally  characterized  with  a 
far-field  test  bed  in  the  laboratory.  Figure  31  shows  the  optical  bench  configuration.  Light  from  an 
argon-ion  laser  at  514  nm  was  directed  through  partially  baffled  optics  to  present  a  collimated  beam 
at  the  retroreflector  array  (RA).  The  returns  from  the  RA  were  reflected  back  through  the  lens 
system,  and  the  array’s  far-field  intensity  distribution  was  recorded  by  a  Pixel  211  digitizing  CCD 
camera  as  shown.  The  Pixel  211  has  12-bit  sampling,  a  192  x  165  format,  and  a  2.46  x  2.46  mm 
active  area.  As  can  be  seen  in  Fig.  31,  the  Fourier  plane  of  the  array  was  located  3.0  m  away  from 
L2.  This  image  plane  was  then  magnified  2.5  times  onto  the  CCD’s  active  area  by  L3  and  L4, 
which  are  achromatic  lenses  selected  to  reduce  aberrations. 

Three  far-field  intensity  distributions  were  captured  as  digitized  images  with  the  CCD  camera. 
Optimum  exposure  time  was  set  to  5.25  s.  The  output  from  the  laser  was  85  mW  and  was  atten¬ 
uated  using  neutral  density  filters  to  prevent  saturation  of  the  CCD.  The  first  FFDP  data  set  was 
generated  with  a  mask  that  exposed  only  one  retroreflector  in  the  16°  incidence  angle  ring.  The 
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second  and  third  data  sets  were  obtained  for  the  full  RA  instantaneous  FFDP  and  time-averaged 
FFDP  respectively. 


6.2  Calibration 

The  CCD  camera  provided  a  digitized  image  of  the  FFDP.  However,  to  relate  the  distribution 
from  raw  pixel  values  to  angular  position  in  //rad  and  ctlrcs  in  m2,  the  sensor  has  to  be  calibrated 
and  conversion  constants  generated. 

6.2.1  Spatial  Calibration 

For  the  spatial  calibration,  the  optical  layout  shown  in  Fig.  31  was  used  with  the  retroreflector 
array  replaced  by  a  circular  aperture  followed  by  a  flat  mirror.  A  circular  aperture/flat  mirror 
combination  observed  at  normal  incidence  has  the  same  high  symmetry  geometry  as  a  normal 
incidence  single  retroreflector,  and  the  FFDP  follows  the  Airy  function  form  [Eq.  (10)].  The  known 
angular  dependence  of  an  Airy  function  served  as  the  basis  of  the  spatial  calibration.  With  the 
mirror’s  normal  aligned  to  the  L2/L3  optical  axis,  an  Airy  pattern  was  projected  onto  the  CCD. 
Images  of  the  projected  Airy  pattern  were  recorded  with  the  CCD  for  12  apertures  with  radius  from 
7  to  18  mm  in  1-mm  steps.  The  weighted  x2  difference  between  the  observed  intensity  distribution 
on  the  CCD  and  an  Airy  intensity  distribution,  Eq.  (10),  was  then  determined  as  a  function  of  the 
scaling  parameter  £  along  a  radial  slice  across  the  image.  The  scaling  parameter  £  converts  the 
CCD  pixel  number  to  an  angle.  The  minimum  weighted  x2  value  then  determines  the  best  fit  Airy 
pattern  for  a  given  aperture,  and  consequently,  relates  a  pixel’s  number  on  the  CCD  xl  to  angular 
position  from 


£xi  =  r—  sin#  .  (18) 

A 

Based  on  this  procedure,  the  resulting  spatial  scaling  coefficient  S  is  0.156  pvdA/pm,  with  a  5.6% 
standard  deviation.  5  also  has  a  detectable  systematic  variation  with  aperture  size  believed  to  be 
due  the  optical  system’s  point  spread  function,  which  was  measured  to  be  ~  70  pm  in  width. 


6.2.2  Radiometric  Intensity  Calibration 

The  CCD  camera  flat  field  correction  and  intensity  calibration  were  determined  by  imaging 
an  integrating  sphere  at  several  power  levels.  Light  from  the  514-nm  laser  was  directed  into  an 
integrating  sphere.  A  UDT  S370  Optometer  measured  the  power  exiting  a  perpendicular  port.  The 
UDT  was  then  replaced  with  the  CCD  camera,  and  an  image  was  recorded.  To  assure  stability 
of  laser  output,  a  second  UDT  monitored  secondary  scatter  from  one  of  the  mirrors  in  the  optical 
train  leading  to  the  integrating  sphere.  Initially,  the  laser  light  intensity  was  set  so  that  the  pixels 
in  the  CCD  camera  all  read  close  to  their  maximum  value  but  were  not  saturated.  The  light  was 
then  attenuated  in  steps  to  the  UDT’s  detection  range  limit.  Images  and  UDT  readings  of  the 
sphere  were  recorded  for  23  different  attenuations. 

For  each  of  the  CCD  camera’s  31,680  pixels,  a  linear  least-squares  fit  yielded 


Pixel  Value(i,  j)  —  g  [td(i,  j )  +  Tl(i,  j)(UDT  Reading)]  , 


(19)- 
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Fig.  32  -  Experimental  ctlecs  for  a  single  retroreflector  at  a  16°  tilt  from  normal 


where  g  is  the  analog  to  digital  gain  and  id{i,j)  and  7 Z(i,j)  are  the  dark  current  and  responsivity 
of  the  i,j  pixel,  respectively.  Equation  (19)  can  be  inverted  to  convert  the  raw  12-bit  data  into 
images  with  units  of  YV/cm2. 

The  relation  ctlrcs  =  yfA2  holds  on  the  optical  axis.  This  determines  the  coefficient  C  that 
scales  the  calibrated  CCD  output  in  W/cm2  to  orrcs  m2-  The  coefficient  C  has  a  mean  value 
of  0.00056,  with  a  standard  deviation  of  22%  for  six  different  aperture  diameters. 


6.3  Results  and  Comparison  to  Numerical  Computation 

The  constants  derived  from  the  calibration  and  bit  map  conversion  discussed  above  were  used 
to  plot  the  FFDP  distributions  (Figs.  32-34).  In  Fig.  32,  the  recorded  far  field  distribution  of 
the  off-axis  single  retroreflector  shows  the  expected  Airy  pattern  elongated  to  some  extent  by  the 
16°  tilt  (see  Figs.  7  and  9).  Figure  33  is  an  instantaneous  view  of  the  entire  array.  The  FFDP 
is  highly  structured  as  the  result  of  interfering  contributions  from  the  array’s  multiple  retroreflec- 
tors.  Atmospheric  turbulence  and  other  effects  cause  time-averaging,  hence,  a  “smoothing”  of  the 
pattern’s  distribution.  Figure  34  shows  a  more  realistic  FFDP  for  the  retroreflector  in  space  when 
pulsed  at  repetition  rates  lower  than  the  rate  of  change  of  the  atmosphere.  This  figure  shows  the 
time-averaged  FFDP  of  the  array  and  corresponds  to  Eq.  (17).  This  figure  shows  how  the  complex 
interference  structure  “washes  out”  and  an  averaged  overall  Airy-like  distribution  prevails. 

Profile  slices  through  the  experimental  and  numerical  values  of  crLRCs  for  a  single  retroreflector 
are  compared  in  the  top  panel  of  Fig.  35.  The  model  locates  the  first  nulls  at  ±79  jurad.  When 
corrected  for  system  effects,  the  experimentally  measured  separation  between  nulls  is  166  ^rad, 
which  represents  a  4.8%  difference  from  the  computed  value.  The  difference  between  numerically 
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Fig.  33  -  Instantaneous  experimental  <tlrcs  for  the  NRL  LEO  array  at  normal  incidence 
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Fig.  34  -  Time- averaged  experimental  ctlrcs  for  the  NRL  LEO  array  at  normal  incidence  (compare  to  Fig.  13) 
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//rad 

Fig.  35  -  Experimented  and  numerical  values  of  <tlrcs  for  single  retroreflector  (top)  at  a  16°  tilt  and  the  NRL  LEO 
array  (bottom)  at  normal  incidence 


predicted  and  experimentally  measured  ctlrcs  was  also  determined.  The  model  predicts  65078m2 
and  the  experimental  ctlrcs  is  52384m2.  The  experimental  value  at  the  peak  is  less  than  the 
numerical  value  by  21.6%. 

For  the  array,  the  same  trends  in  spread  and  magnitude  of  ctlrcs  occur.  As  seen  in  the  bottom 
panel  of  Fig.  35,  the  time-averaged  array  has  no  clear  minimum,  but  the  width  of  the  experimental 
central  maximum  is  wider.  Also,  the  experimental  central  maximum  magnitude  is  lower  than  in 
the  model.  Furthermore,  the  wings  of  the  experimental  <tlrcs  profile  are  higher  than  the  numerical 
prediction.  The  experimental  maximum  for  the  time-averaged  array  ctlrcs  is  272944  m2  while  the 
model’s  maximum  is  319460  m2.  The  difference  is  15%.  At  42  //rad  off  center,  the  experimental 
o'lrcs  again  exceeds  the  numerical  model.  The  measured  <tlrcs  is  122749  m2  and  the  model’s 
o'lrcs  is  124719  m2,  with  a  difference  of  1.6%. 

In  the  context  of  operational  SLR  this  deviation  from  predicted  peak  values  actually  helps.  That 
is,  the  model  calculation  of  ctlrcs  for  a  retroreflector  does  not  include  wave  front  distortion.  It 
does,  as  noted,  include  the  width  of  the  beveled  edges,  the  depth  that  the  retroreflector  is  recessed 
into  the  array,  and  a  scalar  p  for  reflectivity  losses.  The  effect  of  wave  front  distortions  is  to  diffract 
light  from  the  center  maximum  into  relatively  large  angles.  This  enhanced  diffraction  manifests 
as  a  lower  maximum  with  higher  sidelobes.  Wave  front  distortion  has  been  intentionally  included 
in  the  past  for  larger  retroreflectors  (e.g.,  2.5  cm)  and  is  commonly  termed  “spoiling.”  Due  to  the 
velocity  aberration,  the  region  from  39  to  49  //rad  is  of  particular  engineering  significance.  Here, 
the  model  is  very  close  to  the  experimental  values  (Fig.  36). 
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Fig.  36  -  Time-averaged  array  experimental  and  numerical  values  of  <7lrcs  for  left  wing  (top)  and  right  wing  (bottom) 
of  NRL  LEO  array  at  normal  incidence 


7  SPACE  QUALIFICATION  OF  THE  NRL  LEO  RETROREFLECTOR 
ARRAY 


The  NRL  LEO  retroreflector  array  was  qualified  for  space  flight  with  a  number  of  tests.  The 
qualification  unit  is  identical  to  the  flight  articles  except  that  the  bevels  do  not  satisfy  the  0.004-in. 
specification  on  the  retroreflectors.  Because  slight  variations  in  retroreflector  bevel  size  will  have 
no  effect  on  the  mechanical  properties  of  the  array,  substitution  of  retroreflectors  will  not  alter  the 
flight  qualification.  The  flight-qualified  retroreflectors  were  qualified  separately  by  Zygo. 

Data  taken  from  the  tests  described  below  were  analyzed.  No  thermal-related  problems  were 
detected,  and  no  negative  margins  of  safety  for  strength  or  fatigue  were  found.  A  more  complete 
description  of  the  qualification  results  can  be  found  in  Appendix  C. 

The  unit  weighs  7.82  oz  (221  grams),  including  glass,  mount  assembly,  and  bolts  for  mounting. 


7.1  Random  Vibration  Tests 

To  conduct  the  random  vibration  tests,  the  unit  was  mounted  on  an  electrodynamic  shaker. 
The  axial  and  lateral  axes  were  shaken  independently.  Each  vibration  was  2  min.  in  duration.  The 
axial  level  was  13.92  G,  rms.  The  lateral  level  was  also  13.92  G,  rms. 
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7.2  Thermal  Vacuum  Tests 

The  unit  was  mounted  on  a  thermally  controlled  heat  sink,  and  thermocouples  were  used  to 
monitor  and  control  temperature.  The  vacuum  was  less  than  0.00001  torr.  The  unit  was  tested 
through  6  cycles  from  -100°  C  to  +71°  C. 


7.3  Pyroshock  Tests 

The  unit  was  mounted  on  an  electrodynamic  shaker  and  again  the  two  axes  where  shocked 
independently  in  axial  and  lateral  dimensions.  The  peak  levels  of  shock  were  6,000  Gs. 


8  CONCLUSIONS 

Based  on  the  analysis  and  experimental  verification  presented  in  this  report,  we  conclude  that 
the  NRL  LEO  retroreflector  array  will  robustly  close  a  link  from  a  standard  1,100-km  circular  orbit 
using  an  SLR  system  comparable  to  NASA’s  TLRS  above  a  20°  elevation  angle  for  daytime  and 
nighttime  ranging.  The  array  will  also  close  a  link  for  a  system  as  modest  as  the  FTLRS  13-cm 
aperture  system  in  Caparia.  In  addition,  when  analyzed  for  the  meter-class  system  to  be  installed 
at  MRC  in  2000,  the  link  will  be  robust  enough  to  support  unambiguous  ranging  at  even  the  lowest 
elevation  angles,  day  and  night,  without  the  aid  of  an  amplifier. 

The  array  itself  comprises  22  1-cm  diameter  retroreflectors  situated  on  a  hemisphere  to  produce 
a  far  field  diffraction  pattern  that  provides  a  laser  radar  cross  section  of  greater  than  104  m2  for 
elevation  angles  above  20°  with  a  108°  FOV.  It  has  a  mass  of  221  grams  and  measures  82  mm  in 
diameter  by  43  mm  in  height. 

Analysis  was  verified  by  bench  experiments  and  showed  that  manufacturing  defects  in  the  wave 
front  quality  of  the  retroreflector  direct  energy  to  a  small  extent  into  sidelobes.  It  was  shown  that 
these  enhanced  sidelobes  help  the  link  by  compensating  for  velocity  aberration. 

Results  were  compared  to  the  predicted  performance  of  a  single  cube.  It  was  shown  that  the 
FOV  would  be  severely  restricted  when  using  one  cube,  hence  orbital  sampling  would  be  significantly 
reduced. 

Predicted  performance  in  the  near-infrared  at  1064  nm  was  also  presented.  Advantages  of  this 
wavelength  include  better  transmission  through  the  atmosphere  and  covertness  due  to  transmission 
in  the  nonvisible  region  of  the  spectrum. 

The  array  was  space-qualified,  and  testing  specifications  are  detailed  in  Appendix  C. 
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ACRONYMS 


APD  avalanche  photodiode 

COM  center  of  mass 

FAA  Federal  Aviation  Administration 

FFDP  far  field  diffraction  pattern 

FOV  field  of  view 

FTLRS  field-transportable  laser  radar  station 

FWHM  full-width  at  half  maximum 

LEO  low  Earth  orbit 

LRCS  laser  radar  cross  section  (or  ctlrcs) 

MOBLAS  mobile  laser  ranging  system 
MRC  Midway  Research  Center 
NRL  Naval  Research  Laboratory 
OPC  optical  phase  center 
PCA  point-of-closest  approach 
SLR  satellite  laser  ranging 

SOR  Starfire  Optical  Range,  Kirtland  AFB,  New  Mexico 
TLRS  transportable  laser  ranging  system 
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Appendix  A 


DIFFRACTION  PATTERN  AND  PASS  GEOMETRY  CALCULATIONS 


A.l  Circular  Retroreflector  Far  Field  Diffraction  Patterns 

A  convenient  method  to  develop  the  analytic  basis  for  the  numerical  determination  of  far  field 
diffraction  patterns  begins  with  the  simplest  case  and  then  adds  refinements.  With  this  approach, 
numerical  implementation  of  each  stage  can  serve  as  a  limiting  test  case  for  the  next  stage.  Our 
starting  case  is  normal  incidence  without  bevel  losses;  we  finish  with  tilted  incidence  with  bevel 
losses.  As  has  been  noted  in  the  body  of  this  report,  introducing  a  tilt  in  the  incidence  direction  is 
a  much  larger  effect  than  are  the  bevel  losses. 


A.  1.1  Normal  Incidence  Without  Bevel  Losses 


Consider  a  circular  retroreflector  with  bevel  loss  regions  directly  overhead.  The  two-dimensional 
Fourier  transform  of  the  complex  reflectance  of  the  aperture,  effectively  the  retroreflector’s  far  field 
diffraction  pattern  A(x,  y)  is  given  by 


d(kx,ky)  =  ff  dxdyA(x,y)e-ik*xe-ik»y  . 

J  J aperture 


(Al) 


Using  the  function  c(ar)  =  \/V2  -  x2,  which  defines  the  top  half  of  the  circular  aperture,  the  double 
integral  can  be  rewritten  as 


d(kx,  ky)  =  r  dxe~ik*x  f  ^  e-*»vdy 

J—r  J —c(x) 

—  r  e~ikxX  *  \e~ikyc(x)  _  e*fc»c(ar)l  dx 

J—r  - iky  L  -I 

__  r  €-ikxxJ^_  r €ikyc(x)  _ 

J-r  iky  L 


-ikyc(x) 


dx 


—  —  f  e  tkxX  sin  kyc(x) 
foy  J—r 

2  rr 

=  —  I  [cos  kxx  —  i  sin  kxx]  sin  kyc(x)dx 

ky  J—r 

—  f  cos  kxx  sin  ( kyy/r 2  —  x2)  dx 

ky  J  —  r  '  ' 

—  f  sin  kxx sin  ( ky\/r 2  —  x2)  dx 
ky  J  —  r  '  ' 


(A2) 

(A3) 

(A4) 

(A5) 

(A6) 


(A7) 
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When  Eq.  (A7)  is  integrated  numerically,  the  absolute  magnitude  of  a(kx,  ky)  determined,  and  then 
converted  to  a  cross  section,  the  result  is  the  azimuthally  symmetric  Airy  function  of  Fig.  6.  (See 
Section  B.l  for  the  numerical  implementation  of  Eq.  (A7).)  The  numerical  routine  that  implements 
this  geometry  serves  as  a  limit  test  case  for  the  tilted  and  normal  incidence  with  bevel  losses  routines. 


A.  1.2  Tilted  Incidence  Without  Bevel  Losses 


When  a  circular  retroreflector  is  tilted,  the  aperture  changes  shape  to  the  overlap  region  of  two 
offset  circles.  The  tilt  away  from  normal  incidence  can  be  related  to  the  circles’  offset  xq  from  the 
center  of  the  retroreflector.  For  a  retroreflector  of  depth  d  with  a  face  recessed  by  distance  c  and 
index  of  refraction  n,  the  expression  relating  xq  to  the  angle  6inci  between  the  surface  normal  and 
the  incidence  direction  is 

x0  =  c  sin  0inci  +  d  sin  6n  ,  ( A8) 

where 

„  .  /  sin  Bind  \ 

6n  =  arcsin  f — - — J  .  (A9) 

For  a  circular  retroreflector  of  radius  r,  d  =  a/2 r.  The  retroreflector’s  field  of  view  (FOV)  is 
limited  by  the  requirement  that  xq  <  r.  For  a  1-cm  diameter  retroreflector  with  n  =  1.46  and 
c  =  0.1016  cm  (or  0.040  in.),  the  FOV  cut  off  is  at  Qinc{  =  58.579°.  The  aperture  boundary  for  the 
tilted  retroreflector  is  defined  by  c(x)  =  ±y/W~- [x  -  x0)2  for  x  <  0  and  c(x)  =  ±^-  (x  +  x0)2 
for  x  >  0.  The  aperture’s  shape  is  that  of  a  cat’s  iris  or  an  American  football  profile  shown  in  the 
upper  panel  of  Fig.  8. 

Integrating  each  side  of  the  aperture  separately  yields  a  two-dimensional  Fourier  transform  for 
tilted  incidence  with  four  integrals  given  by 


foy)  — 


Y  ^  ^  ^  cos  kxx  sin  [ky\Jr2  -  (x  -  z0)2)  dx 

r(r-x0)  /  - - \ 

+  cos  kxx  sin  (  ky\J r2  -  (x  +  x0)2J  dx 

~  IT  [/  (  ) sin  ^ xX  sin  {^y^r2  ~  (x  ~  xo)2)  dx 

r(r-x0)  /  - - \ 

+  j r  sin  kxx  sin  (  kyyj r2  -  {x  +  x0)2J  dx 


(A10) 


Equation  (A10)  yields  a  reasonable  transform  and  limits  to  an  Airy  function  as  x0  -*  0.  (See 
Section  B.2  for  the  numerical  implementation  of  Eq.  (A10).)  Figure  10  shows  the  sizable  variation 
in  the  laser  radar  cross  section  of  a  tilted  circular  retroreflector  as  one  sweeps  kx  and  ky  with  a 
constant  magnitude  of  k  of  50  fi rad  (i.e.,  roughly  a  satellite’s  velocity  aberration  for  a  1,100-km 
high  orbit). 


A. 1.3  Normal  Incidence  With  Bevel  Losses 

Consider  a  circular  retroreflector  with  bevel  loss  regions  directly  overhead.  There  are  six  regions, 
symmetric  about  the  avaxis,  that  can  be  used  to  compute  the  two-dimensional  Fourier  transform 
for  the  aperture.  These  regions  are  shown  in  Fig.  Al. 
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Fig.  A1  -  On-axis  retroreflector  with  bevel  geometry  (the  six  regions  are  numbered  left  to  right;  regions  3  and  4  are 
each  made  of  two  pieces) 


The  two-dimensional  Fourier  transform  of  the  complex  reflectance  of  the  aperture  A{x,y)  is 
now  given  by 


^ y )  — 


where 


and 


/~x2  rc{x)  f-W 

dxe~lkxX  /  e~lkyVdy  +  /  dxe~lkxX  /  e~tkyVdy 

-r  J—c(x)  J—x  2  Jl+(x) 

r-w/2  [  W”(x)  rc(x) 

+  /  dxe~xkxX  /  e~lkyVdy  +  /  e~lk»ydy 
J—x  1  I  J—c(x)  J—l  (x) 

+  fXl  dxe~ikxX  f~‘  {X)  e~ikyVdy  +  / 

Jw/2  J-c(x)  Jl~ 

L 


c(x) 

l+(x) 


e~%kyV  dy 


rx 2  fl  (*)  fr  rc(x) 

+  I  dxe~tkxX  /  e~'kyVdy  +  /  / 

Jx 2  J—c(x) 


e~tkyVdy 


l+{x)  =  -j=(x  +  w) 


1  ^  =  ' 


(All) 

(A12) 

(A13) 


The  two  integration  limits  along  x  determined  by  the  intersection  of  /+(x)  and  l  (x)  and  the  circular 
aperture  are  given  by  zi  =  \  (—  w  +  \/l2 r2  —  3 w2)  and  X2  =  ^  (ic  +  Vl2r2  —  3 u?2),  respectively. 
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Rather  than  working  with  all  six  terms  in  Eq.  (All)  at  once,  the  forms  are  reduced  to  be 
convenient  for  numerical  integration  individually.  The  six  terms  are: 


(A14) 

(A15) 

(A16) 

(A17) 

(A18) 


^1  ky)  — 

J  X2  dxe~ikxX 

Mx) 

/  e^Hy  , 

J — c(x) 

^2  ( kX  5  ky)  = 

r—w 

/  dxe~ikxX 
J—x  2 

/-*+(*)  .. 

L  *  ■ 

^3  {kxi  ky)  ~ 

r-w/2 

/  dxe~tkx 

J—X\ 

r  [  f 1  (X)  e~ikyVdy  + 
[J-c(x) 

rc(x) 

/-!-(*) 

^4(^X7  ky)  — 

P  dxe~ikxX 
Jw/2 

r~l+(x ) 

/  e~tk*ydy  + 

v/  — c(x)  J 

rc(x) 

/  e 

'»+(*) 

^5(^X7  ky)  — 

j>^l 

^  -iky, 

e  kyVdy  , 

-*-(*) 

e~ikyVdy 


and 


a6(kx,  ky)  =  f  dxe  tkxX  f  ^  J  e~ikyVdy  . 

JX2  J—c(x) 


(A19) 


Both  di(kx,ky)  and  a6(kx,ky)  have  upper  and  lower  limits  defined  by  the  circular  aperture,  as 
in  Eq.  (A2)  for  the  integration  along  y.  Hence,  these  two  terms  can  be  immediately  simplified  to 
the  same  form  as  Eq.  (A7),  yielding 


2  f~x?  /  _ \ 

al(kx,ky)  =  —  J  cos  kxx  sm  [ky\f r2  -  x2)  dx 

2 i  f~x 2  /  , - - 

~  J  sin  kxx sin  (kyVr2  -  x2j  dx 


(A20) 


and 


2  rr  f  _ 

ae(kX)  ky)  =  —  cos  kxx  sin  iky\/r2  -  x 2)  dx 

Ky  Jx 2  '  ’ 

~jrJ  s‘n  kxx  sin  (ky  Vr2  -  x2j  dx 


(A21) 


Note,  of  course,  that  the  x  integration  limits  in  Eqs.  (A20)  and  (A21)  are  different  than  those  in 
Eq.  (A7).  Similarly,  a,2(kx,ky)  and  a^{kXl  ky)  can  be  reduced  to  integrals  over  x  of  the  form 


and 


2  f~w  f  1  \ 

d2(kx,ky)  =  J  cos kxx sin  ycy-j=(x  +  w)J  dx 

2 i  f~w  (  1  \ 

+  J  sin  kxx  sin  ^ky-^=(x  +  w)  J  dx  , 

a5(kx,ky)  —  —  J  cos kxx sin  \^ky—^=(x  -  w) J  dx 

2  i  fx2  (  l  \ 

—  J  sin  kxx  sin  ^=(x  —  w)  J  dx  . 


(A22) 


(A23) 
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Note  the  sign  change  of  the  real  and  imaginary  parts  of  Eq.  (A22). 

Simplifying  ds(kx,ky)  and  a^(kx,ky)  can  be  done  with  the  same  approach.  For  a^{kx.  ky) . 

&^{kxi  ky*)  ” 


ran  T  r— ,+  (x)  •;  fcW  -u 

/  dxe~lkxX  /  e~tk*ydy  +  /  e~tk»ydy 
Jw/2  J—c(x)  Jl+(x) 

_  fXl  c-ikxx  1  _  \eikyl+  (x)  __  eiA:3/c(^)  +  e~^yc(x)  _  e~ikyl+(x)  1  ^ 

*/u>/2  —iky 

=  fXl  e~ikxx  \eikyc(x)  _  e~ikyc(x)  _  c*fcy/+(ar)  _|_  e-*V+(x)J 

Jw  /  2  ky  L 

=  —  /  [*2i  sin  kyc(x)  -  *2i  sin  fej,Z+  (a;) 

^2/  Jw/2  L  ^  ^ 

2  r*1 

—  —  [cos  fcxa;  —  i  sin  fcxa:]  sin  kyc(x)dx 

ky  Jw/ 2 

~  7~  [  [cos ~  i sin A^a:]  sin  (a;)dar 
^2/  Jw/2 

=  i  /  cos  sin  ( ky\J  v2  —  a:2)  dx 
ky  [Jw/2  V  ' 

—  /  cos  fcra:  sin  ( ky-^=(x  +  dx 
Jw/2 


'Vzy 


—  7~  i  f  sin  A:*®  sin  ( ky\fr 2  —  x2\  dx 

ky  [Jw/2  V  J 

—  J  sin  kxxsin  ^ky-^=(x  +  w)^  . 


Finally,  repeating  these  steps  for  63 (kx,  ky)  yields  the  similar  form 

r-w/2 

-an 


/— w/2  f  rc(x) 

dxe~ik*x  /  e~tkyVdy  +  /  e~tkyydy 

-Xi  J—c(x)  J—l  (a:) 


=  t~  I  f  cos  kxx  sin  (kyV r2  —  x2)  dx 
ky  [-/-an  '  ' 

+  j  cosfcxXsin  ^ky-^=(x  —  dx| 

—  7 f  sin  kxx  sin  ( ky\ /  r2  —  a;2]  dx 

ky  (y-an  ' 

+  J  sinfcxxsin  ^ky-^=(x  —  w)^j  dzl 


(A24) 

(A25) 

(A26) 

(A27) 


(A28) 


(A29) 


(A30) 


(A31) 


The  integrals  in  Eqs.  (A20),  (A21),  (A22),  (A23),  (A29),  and  (A31)  are  to  be  evaluated  numer¬ 
ically.  Figure  12  shows  the  variation  in  the  cr lrcs  of  a  beveled  circular  retroreflector  as  one  sweeps 
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kx  and  ky  with  a  constant  magnitude  of  k  of  50  (i rad  (i.e.,  roughly  a  satellite’s  velocity  aberration). 
In  comparison  to  the  variation  seen  in  Fig.  10,  the  bevels  do  not  cause  that  much  change  in  the 
magnitude,  but  they  do  introduce  a  six-fold  symmetry. 


A.  1.4  Tilted  Incidence  With  Bevel  Losses 


For  a  circular  retroreflector  with  bevel  loss  regions,  off-normal  incidence  angles  tilted  about 
the  y- axis  result  in  an  aperture  that  is  a  combination  of  the  cat’s  eye  boundary  and  bevel  loss 
regions.  The  position  of  the  bevel  loss  regions  is  unchanged,  and  Eq.  (A8)  gives  the  correct  offset 
in  this  case  as  well.  Note  that  the  FOV  limit  now  becomes  £0  <  r  —  y.  For  a  1-cm  diameter 
retroreflector  with  n  =  1.46,  c  =  0.1016  cm  (or  0.040  in.),  and  w  =  0.01016  cm  (or  0.004  in.),  the 
FOV  cut  off  is  at  6inc{  =  57.63865°.  Because  the  normal  incidence  case  with  bevel  losses  required 
distinct  integration  regions  on  each  side  of  the  aperture,  the  same  six  regions  used  for  Eq.  (All) 
can  also  be  used  in  this  case.  The  two  integration  limits  along  x  determined  by  the  intersection 
of  ;+(x)  and  l~(x )  and  the  aperture  change  to  xx  =  §  j-  +  *0)  +  ^/|r2  -  |(w  -  z0)2}  and 

x2  —  f  { ( | w  —  £0)  +  7§r2  —  5(10  +  £o)2))  respectively.  The  two-dimensional  Fourier  transform 
for  tilted  incidence  can  be  expressed  in  an  analogous  form  to  Eq.  (All),  where  the  resulting  six 
terms  are  given  by: 


di  {kxi  ky)  — 


Q'<l{kx'>ky)  - 


^3  [kxi  ky  )  — 


^4(^x5  fcy)  — 


J  ^  ^  cos  kxx  sin  (ky^Jr2  -  {x  —  a;0)2)  dx 

~  JT  J  ^  ^  s*n  kx%  sin  -  (x  -  a;0)2)  dx 

-  —  J  cos&^xsin  yky-^=(x  +  w)j  dx 

2 i  [~w  (  1  \ 

+  J  sin  kxx sin  yky—^=(x  +  w)J  dx  , 

jr- 1 J  cos  kxx  sin  [ky^ r2  -  (x  -  x0)2)  dx 

+  J  cosAr^zsin  -  w))  chc| 

2 i  f  f~w/ 2  (  / - \ 

—  t—  \  I  sin  kxx sin  ky\  r2  —  (ar  -  xn)2  dx 

ky  [J-Xl  V  V  J 

r-w/2  /  1  \ 

J  sin  kxx  sin  \  ky~ J=ix  ~  w)j  dx 

[j  f  cos  kxx  sin  (ky\J r2  —  {x  +  x0)2)  dx 

—  J  ^  cos  kxx  sin  ^ky-^=(x  +  w))  dx^ 

-  {y  /  s*n  * xX  s*n  f2  “  (*  +  xq)2)  dx 


+ 

2_ 

k7, 


(A32) 


(A33) 


(A34) 
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and 


fay)  — 


—  J  sin kxxsm  ^kyL=^x  +  w)j  dx|  , 

—  j  cos kxx sin  yky-y^(x  —  w) J  dx 

2 i  rx 2  /  l  \ 

—  —  j  sinkxxsin  yky~^(x  -  w) J  dx  , 


(A35) 


(A36) 


2  f 

ae(kx,ky )  =  —  j  cos 

^2/  «/  a?2 

2i  /•(J—a;o) 


^  ^  &xa:  sin  r2  -  (a;  +  zo)2^ 

2i  r(r-a?o)  f  f - ~\ 

—  /  sin  &x:r  sin  (  kyyr2  —  (a:  +  £o)2  j  cfo 


(A37) 


Equations  (A32)  to  (A37)  limit  properly  to  the  untilted  case  as  zo  — i >  0.  See  Section  B.3  for  the 
numerical  implementation  of  Eqs.  (A32)  to  (A37). 


A. 2  Converting  Pass  Geometry  to  kxkv  Space 

At  a  given  instant,  the  orientation  and  velocity  of  the  satellite  with  respect  to  the  ground  station 
determines  the  effective  aperture  of  the  retroreflector  and  the  velocity  aberration  respectively. 

For  simplicity,  consider  a  nadir- pointing  retroreflector  and  use  a  fixed-Earth  coordinate  system. 
Hence,  the  ground  station’s  location  vector  is  a  constant  denoted  by  G,  with  fixed  spherical  coor¬ 
dinates  of  G,  6gs,  and  <f>gs.  Note  that  0gs  is  the  ground  station’s  colatitude  measured  down  from 
the  north  pole  rather  than  the  latitude  from  the  equator.  The  satellite’s  position  is  denoted  by 
S,  with  spherical  coordinates  of  S,  6S/C,  and  <j>s/c.  The  satellite  velocity  v  includes  the  effect  of 
ground  station  motion.  Finally,  the  range  vector  between  satellite  and  ground  station  is  R  =  S  —  G. 
Figure  A2  illustrates  the  geometry  and  the  local  coordinate  system  at  the  satellite. 

The  angle  of  incidence  which  defines  the  aperture  through  Eq.  (A8),  is  the  angle  between 
S  and  R.  Since  G,  S,  and  R  form  a  triangle,  0snc;  can  be  related  to  the  angular  coordinates  of  G 
and  S  and  the  elevation  angle  0e[ev  of  the  satellite  at  the  ground  site.  The  angle  between  G  and  S 
is  given  by 

cos  7  =  cos  8S/C  cos  6gs  +  sin  9S/C  sin  dgs  cos {<j>s/c  -  4>gs)  .  (A38) 

Hence, 

Oelev  +  2  +  7  +  Qinci  =  ^  (A39) 

or  t r 

Oelev  +  7  +  hnci  =  £  •  (A40) 

In  practice,  the  values  of  0,nct-  and  0eiev  as  a  function  of  time  during  the  pass  are  computed  directly 
by  orbital  dynamics  software.  When  working  with  a  retroreflector  array,  8lrlc.  is  the  angle  between 
R  and  the  array’s  axis.  Hence,  one  must  do  a  little  further  computation  to  determine  the  angle  of 
incidence  for  each  retroreflector  with  respect  to  R. 

The  apparent  velocity  of  the  retroreflector  can  be  expressed  in  terms  of  velocity,  ground  station 
position,  and  satellite  position.  Once  the  velocity  components  in  a  satellite  local  coordinate  system 
are  determined,  the  values  of  kx  and  ky  that  define  the  ground  stations  position  in  the  far  field 
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Fig.  A2  -  Left  panel  shows  the  global  definitions  of  G,  S,  and  R;  right  panel  show  the  coordinate  system  local  to  the 
satellite 


diffraction  pattern  can 

be  obtained  immediately.  The  local  coordinate  system  with 

x  in  the  plane 

defined  by  S  and  R  is 

Z  =  -S  , 

(A41) 

.  RxS 

y  ~  |R  x  S|  ’ 

(A42) 

and 

x  =  y  x  z  . 

(A43) 

Hence 

RxS  _ 

Vy  =  — - -  •  V  , 

|R  x  S| 

(A44) 

(S-G)xS  ? 

|(S  —  G)  x  S|  ' V  ’ 

(A45) 

S  x  G  _ 

|SxG| 

(A46) 

Similarly 

RxS  A  _ 

vx  =  — - ~x—  S-v  , 

|R  x  S| 

(A47) 

[(S-G)xSl 

=  - - - _Lx__s*v  , 

|(S  —  G)  x  S| 

(A48) 

S  x  S  x  G  _ 

|S  x  G| 

(A49) 
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Finally, 

(A50) 

and 

■ 

(A51) 

The  numerical  routine  to  determine  kx  and  ky  from  the  results  provided  by  the  orbital  dynamics 
group  is  a  distinct  module  to  simplify  testing.  (See  Section  B.4  for  the  numerical  routine  that 
determines  kx  and  ky  from  the  results  provided  by  the  orbital  dynamics  group.) 


Appendix  B 


SCILAB  LISTINGS 


B.l  Circular  Retroreflector  Directly  Overhead  Without  Bevels 

function  [  z  ]  =  curve (  theda.phi  ) 

erra=10“(-16) 

errr=10“(-12) 

theda=theda*10“ (-6) 

r=0.005 

kl  =  2*7.pi/(0. 532*10“  (-6)) 
kx=kl*cos (phi) *theda 
ky=kl*sin(phi) *theda 

Ireal=2.0*integrate( [’if  abs(ky)  >  (10“(-6))  then  cos(kx*x)* 
sin(ky*sqrt(r“2-x“2))/ky, ’ ; ’else 
cos(kx*x)*sqrt(r“2-x“2) .end’] , ’x* ,-r,r,erra,errr) 
Iimg=2.0*integrate([’if  abs(ky)  >  (10“(-6))  then  sin(kx*x)* 
sin(ky*sqrt (r“2-x“2) ) /ky , ’ ; ’ else 
sin(kx*x)*sqrt(r“2-x“2) .end’] , ’x’ ,-r,r,erra,errr) 
z=4.0*,/,pi*(Ireal“2+Iimg“2)/(  (0.532*10“ (-6)) “2  ) 

B.2  Circular  Retroreflector  Off  Normal  Incidence  Without  Bevels 

function  [  z  ]  =  lr (alpha, theda, phi  ) 

erra=10“(-16) 

errr=10“(-12) 

theda=theda*10“ (-6) 

n=l .46 

c=0 .0001 

alphan=as in ( s in ( alpha) /n) 

r=0 . 005 

r2=r“2 

d=2*r*sqrt(2) 

off set=c*tan (alpha) +d*tan(alphan) 
kl  =  2*’/,pi/(0. 532*10“  (-6)) 
kx=kl*cos (phi) *theda 
ky=kl*s in (phi ) *theda 

llreal=2.0*integrate( [’if  abs(ky)  >  (10“(-6))  then  cos(kx*x) 
*sin(ky*sqrt(r2-(x-offset)"2))/ky, ’ ; ’else 
cos(kx*x)*sqrt(r2-(x-offset)“2) .end’] , ’x’ ,-(r-offset) , 
0.0,erra,errr) 
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ulreal=2.0*integrate([,if  abs(ky)  >  (10“(-6))  then  cos(kx*x) 
*sin(ky*sqrt(r2-(x+offset)“2))/ky, * ; ’else 
cos(kx*x)*sqrt(r2-(x+off set)“2)  ,end’]  ,  ’x’  ,0.0,  (r-offset)  , 
erra,errr) 

llimg=2.0*integrate( [’ if  abs(ky)  >  ( 10~ (—6) )  then  sin(kx*x) 
*sin(ky*sqrt(r2-(x-offset)“2))/ky, ’ ; ’else 
sin(kx*x)*sqrt(r2-(x-off set)“2)  ,end’]  ,  ’x’ , -(r-offset) , 
0.0,erra,errr) 

ulimg=2 .0*integrate(  [’ if  abs(ky)  >  (10“(-6))  then  sin(kx*x) 
*sin(ky*sqrt(r2-(x+offset)“2))/ky, ’ ; ’else 
sin(kx*x)*sqrt(r2-(x+off set)“2) .end’] , ’x’ ,0.0, (r-offset) , 
erra.errr) 

z=4 . 0*‘/,pi*  ( (llreal+ulreal)  “2+  (llimg+ulimg)  “2)  /  (  (0 . 532*  10"  (-6) )  "2  ) 


B.3  Circular  Retroreflector  Off  Normal  Incidence  With  Bevels 
function  [  Ires  ]  =  kLRCS(  theda,K,k) 

//  calculates  the  laser  radar  cross  section  for  a  retro,  corner  cube. 

//  The  cube  has  a  round  aperture  of  radius  r  and  beveled  edges  of 
//  width  w.  It  is  recessed  a  depth  c  in  the  holder.  Its  index  is  n. 

//  INPUT:  theda  in  radians,  k  and  K  also. 

//  theda  is  the  tilt  angle  and  K  and  k  are  the  x  and  y  angles  off  the 
//  incident  angle. 

//  This  avoids  k=0  singularity,  k  implicitly  cancels  if  less  than  10“(-6). 
if  abs(k)  <  (10“(-14))  then  k=(10“ (-14) ) ,  end 

ka=l . 1810498*10~7  //  2pi/532.0*10"(-9) ; 

k=ka*k 

K=ka*K 

ro=3.33*10"(13)  //  0.75*4*pi/(0.532*10“(-6))"2 
r=0.005  ;  //  cube  aperture  radius. 

c=0.0015  //  1.0*10“ (-3)  ;  //  depth  of  recess  10“-3  is  about  40  mils 

d=0. 00707  ;  //  distance  to  vertex,  a*sqrt(2) 

w=(1.016e-4)*2  //  width  of  camfer  on  the  cube’s  corner,  4  mills=l .016e-4 

n=1.46  ;  //  index  of  refraction 

thedan=asin(sin(theda)/n) ; 
offset=c*tan(theda)+d*tan(thedan) ; 
a=tan('/,pi/6) 
b=w/  (2*cos  (*/,pi/6) ) 
xl=w/2; 

x2=w  //  w/(2*sin(30degs))=w; 
alp=a“2+l 

bet=2*a*b+2*off set 
gam=b“2+offset“2-r“2 

x3=real (-bet+sqrt (bet“2-4*alp*gam) ) / (2*alp) 
bet=-2*a*b+2*off set 
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x4=real (-bet+sqrt (bet ~2-4*alp*gam) ) / (2*alp) 
x5=r-off set 

if  (  x2  <  x4  )  &  (  xl  <  x5  )  fe  (  x4  <  x5  )  &  (  xl  <  x3  )  then 
kr=k*r; 

zl=acos((xl+offset)/r) ; 
z3=acos((x3+offset)/r) ; 
z4=acos ( (x4+of f set) /r) ; 
z5=real(acos((x5+offset)/r)) ; 
y4= (x4+of f set ) /r ; 
y5=1.0  //  (x5+off set)/r; 
mzl=acos((-xl-offset)/r) ; 
mz3=acos((-x3-offset)/r) ; 
mz4=acos((-x4-offset)/r) ; 
mz5=real(acos((-x5-offset)/r)) ; 
ak=a*k; 
bk=b*k; 

//  Substitution:  cos(z)=(x+offset)/r  for  even  integrals  above  II.  This 
//  helps  the  argument  of  the  square  root. 

//  It  is  written  for  better  convergence. 
errl=10~ (-12) ; 
err2=10~ (-8) ; 

//  upper  boundary  of  II 

Il=integrate(  ,sin(ak*x-bk)*cos(K*x)/k’  , ’x’,x2  ,  x4  ,errl,err2  ), 

//  upper  boundary  of  III 

I2=integrate(  ’-r+sin^r+sqrtQ.-y'^^+cosOC’Kr+y-offset))/^  ,  *y*  ,y4,y5, 
errl,err2  ), 

//  lower  boundary  of  I 

I3=-integrate(  ’sin(ak*x+bk)*cos(K*x)/k’  , ’x’ ,xi,x3,errl,err2  ), 

//  upper  boundary  of  I 

I4=integrate(  ,r*sin(z)*sin(kr*sin(z))*cos(K=^=(r*cos(z)-offset))/k,  ,  >z> , 
zl,  z3,errl,err2  ), 

//  upper  boundary  of  II 

I5=integrate(  ’sin(ak*x-bk)*sin(K*x)/k’  ,,x,,x2  ,  x4  ,errl,err2  ), 

//  upper  boundary  of  III 

I6=integrate(  ’r*sin(z)*sin(kr*sin(z))*sin(K*(r*cos(z)-offset))/k’  ,  >z’ , 
z4,z5  ,errl,err2  ), 

//  lower  boundary  of  I 

I7=-integrate(  ,sin(ak*x+bk)*sin(K*x)/k>  ,’x’,  xl,  x3,errl,err2  ), 

//  upper  boundary  of  I 

I8=integrate(  ,r*sin(z)*sin(kr*sin(z))*sin(K*(r*cos(z)-offset))/k>  , ’z’ , 
zl,  z3,errl,err2  ), 
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//  upper  boundary  of  -II  ;  -II  is  the  left  mirror  image  of  II 
I9=integrate(  J sin(-ak*x-bk)*cos(K*x)/k>  , ’x’,^  ,-x2  ,errl,err2  ), 

//  upper  boundary  of  -III 

I10=integrate(  >r*sin(z)*sin(kr*sin(z))*cos(K*(r*cos(z)+off set))/k>  , ’z> , 
mz5 , mz4 , err 1 , er r2  ) , 

//  lower  boundary  of  -I 

Ill=-integrate(  ’sinC-ak+x+bk^cosCK+x)/^  , >x’,-x3  ,  -xl  ,errl,err2  ), 

//  upper  boundary  of  -I 

I12=integrate(  ’r+sinCz^sinCkr+sinCz^+cosCK+Cr+cosCzHoffset))/^  , 'z* , 
mz3,  mzl,errl,err2  ), 

//  upper  boundary  of  -II 

I13=integrate(  ’sinC-ak+x-bk^sinCK+x)/^  ,***  9-xl  ,  -x2  ,errl,err2  ), 

//  upper  boundary  of  -III 

I14=integrate(  ’r+sinCz^sinCkr+sinCz^+sinCK+Cr+cosCz^offset))/^  ,  ,z> , 
mz5,mz4,errl,err2  ) , 

//  lower  boundary  of  -I 

I15=-integrate(  * sin(-ak*x+bk)*sin(K*x)/k*  ,  JxJ,  -x3,  -xl,errl,err2  ), 

//  upper  boundary  of  -I 

I16=integrate(  ’r+sinCz^sinCkr+sinCz^+sinCK+Cr+cosCz^offset))/^  , , 
mz3,mzl,errl,err2  ) , 

//  combine  all  contributions 

Ires  =  (2)*  (1 1-I2+I3-I4+I9-1 10+111-112-*/, i*  (I5-I6+I7-I8+1 13-1 14+1 15-1 16)  ), 

else  lrcs=0,end 

lrcs=ro*abs (lrcs*conj (Ires) ) 


B.4  Conversion  to  kxky 


mode(-l) 

//  Program  generates  LRCS  for  passes  found  in  sub  directory  passes. 
// 

//  cube  init,  enum  is  the  number  of  cubes  to  be  considered. 

unix_w( ’date* ) 

cnum=22; 

print  (*/,io(2)  ,cnum) 

//  thada  is  the  sperical  angle  from  the  z  axis. 

//  phi  is  the  angle  from  the  x  axis. 

c=2 . 99*10~8 

thada=2**/,pi; 


cub=zeros(cnum,3) ; 
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phi=16*'/,  pi/180  ; 

for  i=l:4,cub(i, l)=cos( (i-l)*thada/4)*sin(phi) ;end 
for  i=l :4,cub(i,2)=sin((i-l)*thada/4)*sin(phi) ;end 
for  i=l:4,cub(i,3)=cos(phi) ;end 
phi=32*'/,pi/ 180 ; 

for  i=8: 12,  cub (i,  l)=cos(*/,pi/8+(i“l)*thada/8)*sin(phi)  ;end 
for  i=5: 12,cub(i,2)=sin(#/lpi/8+(i”l)*thada/8)*sin(phi)  ;end 
for  i=5 : 12 , cub(i,3)=cos(phi) ; end 
phi=48**/.pi/ 180 ; 

for  i=13 : 22, cub(i, l)=cos((i-l)*thada/10)*sin(phi) ;end 
for  i=13:22,cub(i,2)=sin((i-l)*thada/10)*sin(phi) ;end 
for  i=13:22,cub(i,3)=cos(phi) ;end 

//  klrcs  returns  LRCS  for  one  cube  kLRCS(theda,kxc ,kyc) ,532nm 
; getf ( " . /klrcs . sciM , Mc" ) ; 

//  Rot  generates  the  rotation  matrix  to  move  the  array  coordinates 
//to  the  cube  coordinates. 

//In  function  kLRCS  from  klrcs. sci  the  major  axis  of  the  cube  aperture 
//  is  always  the  y  axis. 

//  Rxc  normalized  is  the  major  axis  and  is  pi/2+phi  ,  phi  the  rotation  axis. 
//  cos(pi/2+phi)=-sin(phi) 
deff ( ’ [Rc]=Rot (R,c) ’ , [ 

,Rcr= [R(2)*c(3)-R(3)*c(2);R(3)*c(l)-R(l)*c(3);R(l)*c(2)-R(2)*c(i)] » ; 
>n=sqrt(Rcr(l)*Rcr(l)+Rcr(2)*Rcr(2)+Rcr(3)*Rcr(3)) ’ ; 

*if  n  —  0  then  Rc=[l , 0; 0 , 1]  , 

else  Rc=[Rcr(2)/n,-Rcr(l)/n;Rcr(l)/n,Rcr(2)/n] jend* 

]) 

deff ( } [z]=cross (x,y) » , ,z=[x(2)*y(3)“x(3)*y(2) ,x(3)*y(l)-x(l)*y(3) , 
x(l)*y(2)-x(2)*y(l)] >) 

II  macros  to  extract  data  from  pass  geometry  found  in  a  file  in  passes/ 
deff ( *  [r]=r(i) 9  , ,r=y(i,32) }) 

deff ( * [R]=R(i) , , ,R=[y(i,18),y(i,19),y(i,20)] J) 
deff ( 9  [s]=S(i) }, ,s=[y(i,6),y(i,7),y(i,8)]  0 
deff ( 9  [g] =G(i) } , ,g=~ [y(i,12),y(i,13),y(i,14)] ’ ) 

deff ( 7  [v]=V(i) J , ,v=-[y(i,9),y(i,10),y(i,ll)]-[y(i,15),y(i,16),y(i,17)]  J) 
deff ( 9 [ele]=el(i) , , ,ele=y(i,31) >) 

deff ( J [t]=time(i) J , ’t=y(i,3)*3600+y(i,4)*60+y(i,5) ») 

// 

//  LRCSA  returns  the  LRCS  for  the  22  cube  array 

II 

deff ( ’ [lra]=LRCSA(theda,kxa,kya,cub) * , [  } lra=0 . 0  * ; 

,R=[sin(theda) ,0, cos (theda)] * ; 

’for  i=l:cnum,  c=[cub(i, 1) ;cub(i,2) ;cub(i,3)] ;  kc=Rot (R, c)* [kxajkya] ; 
lra=lra+kLRCS(acos(cub(i ,3) ) ,kc(l) ,kc(2)) ;  end’]  ) 

// 

//  get  the  data  for  passes 

// 

II 

II  Begin  of  main  loop,  reading  in  data  from  a  different  file 

// 

for  p=l:l ,//  filnum(l,l).  This  was  for  multiple  passes.  Option  is  abandoned, 
f ilename= ’ lpm03234  * 
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print ('/.io (2) ,  filename) 
load (filename) 

//  open  file  for  output,  t ime, Ires , elevation 
funit=f ile( ’open* , 'array .dat ' , 'new' ) 
ysize=size(y) 

//  creating  arrays  for  storing  data  for  future  use. 
lc=zeros (22 , ys ize (1,1)) 

1=1 :ysize(l,l) 

1=1*0 

kx=l 

ky=l 

t=l 

ele=t 

vel=t 

ang=t 

slantr=t 

for  j=l:ysize(l, 1) , 

//  calculating  the  incident  vector,  normalized,  and  in  the  array  coordinates 
nR=R( j )/norm(R( j ) ) 
nS=S(j )/norm(S( j )) 

//  z  axis  of  the  array  is  parallel  to  S,  nadir  pointing.  Nza=cos(theda_inc) 
Nza=-nR*nS ' 

sintheda=sqrt ( 1-Nza~2) 

//  Nadir  case  has  Nza=l  and  RxS=0,  nRxS=  INF 
if  Nza  <1.0  then 
RxS=cross (R( j ) ,S( j ) ) ; 
nRxS=RxS/norm(RxS) ; 
nV=V ( j ) /norm (V ( j ) ) ; 

Nya=sintheda*(nRxS*nV' ) ; 

RxSxS=cross (RxS , S ( j ) ) ; 
nRxSxS=RxSxS/norm(RxSxS) ; 

Nxa=sintheda* (nRxSxS*nV ' ) ; 

Ninc=[Nxa,Nya,Nza]  ; 
else  Ninc=[0,0, 1] ,end; 

//  V  perpendicular  to  the  direction  of  incidence 
VxR=cross(V(j) ,R(j ) )/norm(R( j ) ) ; 

Vp=norm(VxR) ; 

//  va,  velocity  aberation 

va=2*Vp/c 

vel(j )=va 

//  project  each  cube  on  to  the  plane  perpendicular  to  R  then  calculate  Ires 
//  test  one  nadir  pointing  cube 
for  i=l:cnum, 

Nc=[cub(i, 1) ,cub(i,2) ,cub(i,3)] ; 
theda=acos ( (Ninc*Nc 9 ) ) ; 

//  The  next  is  parallel  to  the  cubes  y  axis 
yc=cross(Ninc,Nc) 

Nyc=yc/norm(yc ) 
vay=va*Nyc(l) 
vax=va*Nyc(2) 
kx(j)=vax 
ky(j)=vay 
ang(j)=theda 

if  theda  <  0.9  then  lc(i,  j )=kLRCS(theda,vax,vay)  ,end 


l(j )=1( j )+lc(i, j ) 

,end 

ele(j)=y(j,31) 

t(j)=time(j) 

slantr ( j )=norm(R( j ) ) 

//  End  of  loop  for  one  pass 
,end 

//  write  time, elevation,  and  Ires  to  a  file 

for  j=l :ysize(l , 1) , write (f unit , [t (j )  ,l(j ) ,ele(j ) , slantr(j )] , * (4el2. 5) 7 ) ,end 
f ile(Mclose" ,funit) 

//  save  Ires  for  each  point/cube 

save(,221rcs> ,1c) 

save( >22t , ,t) 

save( ^ele’ ,ele) 

save( ,22range) ,slantr) 

xsetech ( [0 , 0 , 1 . 0 , 0 . 5] ) 

plot (t , 1 ) 

xtitleC’LRCS  vs  time*) 
xsetech ([0,0.5, 1.0, 0.5]) 
plot(t,ele) 

xtitle( 3 Elevation  vs  time’) 
xsetech( [0 , 0 , 1 . 0 , 1 . 0] ) 
pca=0 . 0 

for  ipca=l:ysize(l,i) , 

if  pea  <  y (ipca, 31)  then  pca=y(ipca, 31) , 

date= 9 September  >+string(y (ipca, 1) )+> , ,+string(y(ipca,2) )+J 
,+string(y(ipca,3) )+’ : >+string(y(ipca,4) ) 

, end , end 

scnum=str ing ( enum) 

title=string(pca)+>  degree  PCA,+>  ,+date+»  c.number  'tsennm 
xtitle(title) 

//  End  of  main  loop 
,end 

unix_w( ’date  * ) 


52 


Appendix  C 


SPACECRAFT  ENGINEERING  DEPARTMENT  SPACE  QUALIFICATION 

REPORT 


C.l  Test  Report  for  the  Retroreflector  Array 

August  1996 
Robert  B.  Patterson 
Design  Section,  Code  8211 
U.S.  Naval  Research  Laboratory 


C.1.1  Random  Vibration  Tests 

The  test  article  (RR-ST-1100-01  assembly)  was  subjected  to  two  random  vibration  tests.  One 
in  the  axial  (Z-axis)  direction  and  one  in  the  lateral  (X-axis)  direction.  Since  the  assembly  is  a 
symmetrical  round  dome,  a  test  in  the  Y-axis  direction  was  considered  unnecessary. 

The  test  article  was  assembled  using  flight-quality  parts  with  the  exception  of  the  optics,  which 
are  identical  to  the  flight-quality  units  in  all  physical  characteristics,  but  have  slight  flaws  in  the 
optical  properties.  The  qualification  dome  (RR-ST-1030-1)  was  not  irradiated;  this  will  have  no 
effect  on  the  fit  of  the  parts  to  one  another  or  the  strength  or  stiffness  of  either  that  part  or  the 
assembly. 

The  lateral  test  was  run  first.  The  accelerometer  data,  from  an  accelerometer  mounted  on  the 
shake  table  close  to  the  test  article,  are  shown  in  Fig.  Cl.  The  test  article  survived  the  random 
vibration  test  in  the  lateral  axis. 

The  axial  test  was  run  next.  The  accelerometer  data  are  shown  in  Fig.  C2.  The  test  article 
survived  the  random  vibration  test  in  the  axial  axis. 


C.l. 2  Shock  Testing 

The  shock  tests  were  next.  Again,  only  an  axial  and  lateral  test  were  needed.  The  test  was 
performed  on  the  same  shake  table  that  was  used  to  perform  the  random  vibration  tests.  The 
accelerometer  data  were  obtained  the  same  way,  with  the  accelerometer  mounted  on  the  test  fixture 
close  to  the  test  article. 

The  axial  shock  test  was  completed  first.  The  test  article  passed.  The  accelerometer  data  are 
shown  in  Fig.  C3. 

The  lateral  shock  test  was  then  run.  The  test  article  passed.  The  accelerometer  data  are  shown 
in  Fig.  C4. 
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Fig.  Cl  -  Lateral  accelerometer  data 
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Fig.  C2  -  Axial  accelerometer  data 
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Fig.  C3  -  Axial  shock  accelerometer  data 
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Fig.  C4  -  Lateral  shock  accelerometer  data 
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TIME  (MINUTES) 

Fig.  C5  -  Thermal  vacuum  test  temperature  profile 


C.1.3  Thermal  Vacuum  Testing 

The  qualification  assembly  was  placed  in  a  small  vacuum  chamber  and  cycled  between  the 
survival  temperature  extremes  to  verify  that  the  assembly  could  survive  the  environments  it  will 
see  while  attached  to  the  host.  This  was  done  by  bolting  the  assembly  to  a  small  “cold”  plate  that 
had  its  temperature  controlled  by  allowing  liquid  nitrogen  to  flow  through  copper  tubing  that  was 
soldered  to  the  backside  of  the  copper  “cold”  plate  or  by  applying  power  to  electrical  heaters  that 
are  attached  to  the  front  side  of  the  plate,  or  some  combination  of  the  two.  The  plate  had  two 
thermocouples  mounted  to  it.  One  was  for  feedback  to  the  controller  used  to  control  the  flow  of 
liquid  nitrogen;  the  other  was  used  as  the  data  for  the  test. 

The  test  article  was  placed  in  the  vacuum  chamber,  which  was  then  pumped  down  to  a  vacuum 
of  at  least  0.0001  torr.  The  assembly  was  then  subjected  to  temperature  extremes  of  71°C  for  the 
high  temperature  and  ~-100°C  for  the  low  temperature.  The  assembly  was  held  at  each  temperature 
extreme  for  a  1-h  dwell  period  to  ensure  that  the  test  article  stabilized  at  the  test  temperature. 
The  test  consisted  of  6  cold  dwells  alternated  with  6  hot  dwells. 

The  qualification  assembly  passed  the  test  without  suffering  any  failures.  The  test  went  accord¬ 
ing  to  plan,  except  when  the  control  valve  stuck  open  and  the  temperature  went  down  to  —  157°C 
during  the  second  cold  dwell.  Since  the  temperature  stayed  below  the  maximum  of  —  100°C,  this 
did  not  negate  this  part  of  the  test.  The  faulty  valve  was  changed  out  during  the  next  transition 
from  cold  to  hot  dwell  and  the  test  continued.  A  time  vs  temperature  graph  of  the  test  profile  is 
shown  in  Fig.  C5. 
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C.2  Test  Procedure  for  the  Retroreflector  Array 

June  19,  1996 
Robert  B.  Patterson 
Design  Section,  Code  8211 
U.S.  Naval  Research  Laboratory 


C.2.1  Test  Definition 

C.2. 1.1  Scope  This  document  contains  the  requirements  and  procedures  for  the  random  vi¬ 
bration,  shock,  and  thermal  testing  to  be  performed  on  the  retroreflector  array  in  order  to  qualify 
it  for  space  flight  on  its  host. 

C.2. 1.2  Test  Objective  The  objective  of  this  testing  is  to  ensure  that  the  retroreflector  array  is 
structurally  sound  and  can  withstand  the  launch  and  postlaunch  environments  without  experiencing 
structural  failure. 


C.2. 1.3  Pass/Fail  Criteria 

1.  No  parts  (spring,  retaining  ring,  retro  seat,  or  retro  optic)  ejected  from  the  assembly. 

2.  No  noticeable  permanent  deformation  around  the  seat  area  at  the  optic  opening  in  the  retrore¬ 
flector  housing. 

3.  No  noticeable  permanent  deformation  around  the  retaining  ring  groove  area  of  the  retrore¬ 
flector  housing. 

C.2. 1.4  Test  Description  Since  the  retroreflector  assembly  is  a  round  part  that  is  symmetrical 
about  the  axis  of  rotation,  we  will  only  have  to  test  it  in  one  longitudinal  direction  and  the  axial 
direction.  We  will  mount  the  retroreflector  assembly  to  the  shaker  table  so  that  it  is  oriented  for 
testing  in  the  longitudinal  direction.  We  will  then  run  the  test  article  through  the  13.92  g  rms 
random  vibration  test  environment.  We  will  then  “hit”  the  test  article  with  an  extremely  sharp 
impulse  to  simulate  the  shock  environment  (6000  g)  specified  by  the  host  vehicle.  We  will  now 
remount  the  retroreflector  assembly  so  that  it  is  oriented  in  the  axial  direction  and  repeat  the 
random  vibration  and  shock  tests. 

The  thermal  test  consists  of  bolting  the  retroreflector  assembly  to  a  “cold”  plate  (that  also 
has  heaters  installed  on  it),  placing  this  assembly  into  a  vacuum  chamber,  and  then  cycling  the 
assembly  through  its  survival  temperature  extremes  six  times. 

C.2. 1.5  Test  Location  The  test  will  be  conducted  in  the  Naval  Center  for  Space  Technology’s 
Payload  Checkout  Facility,  Building  A-59,  at  the  Naval  Research  Laboratory  in  Washington,  DC. 

C.2. 1.6  Responsibilities  The  Spacecraft  Engineering  Department  is  responsible  for  the  over¬ 
all  performance  of  the  test,  the  data  reduction,  and  production  of  the  final  test  report. 
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C.2.1.7  Test  Log  The  test  conductor  shall  maintain  a  daily  test  log  that  shall  include  a  detailed 
record  of  the  test  progress,  significant  events  and  any  deviations  from  the  test  procedures  outlined 
in  this  document. 


C.2.1.8  Photographic  Coverage  Photographs  will  be  taken  of  the  test  specimens,  test  fix¬ 
tures,  and  peripheral  test  equipment.  Each  transducer  location  will  be  marked  and  photographed. 


C.2.2  Test  Conditions 

C.2.2.1  General  The  test  will  be  conducted  in  the  ambient  temperature  and  pressure  in  the 
Payload  Checkout  Facility,  Building  A-59,  Naval  Research  Laboratory. 


C.2.2. 2  Applied  Loads  The  test  articles  shall  be  tested  to  the  levels  described  by  the  host 
vehicle  as  listed  in  Appendix  A  NCST  Qual.  Procedures.  The  qualification  unit  shall  be  tested  to 
qualification  levels,  and  the  flight  unit  shall  be  tested  to  acceptance  levels. 


C.2.2.3  Test  Report  A  final  test  report  will  be  prepared  to  document  the  test  results.  This 
report  shall  include  overall  test  objectives,  test  article  and  test  setup  descriptions,  reduced  test 
data,  test  anomalies  and/or  significant  observations,  and  a  summary  of  the  test  results. 


C.3  Final  Acceptance  Test  Report  of  the  Retroreflector  Array 

November  1996 

Robert  B.  Patterson 

U.S.  Naval  Research  Laboratory 


C.3.1  Random  Vibration  Tests 

The  flight  units  (RR-ST-1 100-01  assembly)  were  each  subjected  to  two  random  vibration  tests, 
one  in  the  radial  (Z-axis)  direction  and  one  in  the  lateral  (X-axis)  direction.  Since  the  assembly  is 
a  symmetrical  round  dome,  a  test  in  the  Y-axis  direction  was  not  considered  necessary. 

The  flight  units  were  assembled  using  flight-quality  parts  and  were  labeled  assembly  A,  B,  and 
C.  The  flight  units  were  tested  to  acceptance  (flight)  levels  6  dB  below  qualification  for  a  total  of 
6.92  g  rms. 

The  lateral  test  was  run  first.  The  accelerometer  data,  from  an  accelerometer  mounted  on  the 
shake  table  close  to  the  test  article,  are  shown  in  Figs.  C6,  C7,  and  C8.  The  flight  units  survived 
the  random  vibration  test  in  the  lateral  axis. 

The  axial  tests  were  run  next.  The  accelerometer  data  are  shown  in  Figs.  C9,  CIO,  and  Cll. 
The  flight  units  survived  the  random  vibration  test  in  the  axial  axis. 
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Fig.  C12  -  Axial  shock  accelerometer  data,  acceptance  level  S/n  -  A 


C.3.2  Shock  Testing 

The  shock  tests  were  next.  Again,  only  an  axial  and  lateral  test  were  needed.  The  test  was 
performed  on  the  same  shake  table  that  was  used  to  perform  the  random  vibration  tests.  The 
accelerometer  data  were  obtained  the  same  way,  with  the  accelerometer  mounted  on  the  test  fixture 
close  to  the  test  article. 

The  axial  shock  test  was  completed  first.  The  flight  units  passed.  The  accelerometer  data  are 
shown  in  Figs.  C12,  C13,  and  C14. 

The  lateral  shock  test  was  then  run.  The  flight  units  passed.  The  accelerometer  data  are  shown 
in  Figs.  C15,  C16,  and  C17. 
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Fig.  C17  -  Lateral  shock  accelerometer  data,  acceptance  level  S/n  -  C 


